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Abstract: MicroRNAs, which are small endogenous RNA regulators, have been associated 
with various types of cancer. Breast cancer is a major health threat for women worldwide. 
Many miRNAs were reported to be associated with the progression and carcinogenesis of 
breast cancer. In this study, we aimed to discover novel breast cancer-related miRNAs 
and to elucidate their functions. First, we identified confident miRNA-target pairs by 
combining data from miRNA target prediction databases and expression profiles of 
miRNA and mRNA. Then, miRNA-regulated protein interaction networks (PINs) were 
constructed with confident pairs and known interaction data in the human protein reference 
database (HPRD). Finally, the fiinctions of miRNA-regulated PINs were elucidated by 
functional enrichment analysis. From the results, we identified some previously reported 
breast cancer-related miRNAs and functions of the PINs, e.g., miR-125b, miR-125a, 
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miR-21, and miR-497. Some novel miRNAs without known association to breast cancer 
were also found, and the putative functions of their PINs were also elucidated. These 
include iniR-139 and iniR-383. Furthermore, we validated our results by receiver operating 
characteristic (ROC) curve analysis using our miRNA expression profile data, gene 
expression-based outcome for breast cancer online (GOBO) survival analysis, and a literature 
search. Our results may provide new insights for research in breast cancer-associated miRNAs. 

Keywords: miRNA; breast cancer; protein interaction network; functional analysis 



1. Introduction 

Breast cancer is a global health threat for women. According to a 2008 survey [1], breast cancer 
was the leading cause of cancer deaths in women. Our knowledge of possible risk factors has led to 
developments in diagnostic methods, drugs, and surgery procedures for treatment [2,3]; however, the 
details of breast carcinoma progression, and perhaps most importantly, how to cure breast cancer, 
remain elusive. 

Previous research has identified a number of risk factors for breast cancer. Early menarche, late 
menopause, obesity, late first full pregnancy, and hormone replacement therapy were considered as 
high risk factors for breast cancer [2]. Breast cancer risk has also been reported to be related to fat 
intake in diets rich in red meats and high-fat dairy foods [3]. 

MicroRNAs (miRNAs) [4], are short endogenous non-coding RNAs which are able to regulate gene 
expression. After miRNA precursors are transcribed fi-om the genome or generated fi-om spliceosomes, 
they are exported to the cj'toplasm and further processed by the Dicer complex [5]. The mature miRNA is 
then bound to Argonaute protein, forming a miRNA-protein complex known as the RNA-induced 
silencing complex (RISC) [6], miRNP, or RNAi (RNA interference) enzyme complex [7,8]. The RISC 
has been reported to down-regulate target genes by translational repression [9] or mRNA cleavage [10]. 

Like other protein-based regulators, miRNAs have been associated with cancer. Calin et ah, 
reported that miR-15 and miR-16 were deleted in leukemia [11], which was believed to be one of the 
earliest reports associating miRNAs with cancer [12]. After this report, many miRNAs were found to 
act as tumor suppressors or oncogenes (also known as oncomirs). For example, miR-21 was identified 
as an oncomir in hepatocellular cancer [13], breast cancer [14], and kidney cancer [15]. On the other 
hand, let-7c was found to be a tumor suppressor in prostate cancer [16], and miR-181a was reported as 
a tumor suppressor in glioma [17]. Further, miR-125b [18], and miR-145 [19] were identified as tumor 
suppressors in breast cancer, and miR-125a was found to repress tumor growth in breast cancer [20]. 
Thus, it is highly likely that miRNAs play an important role in breast cancer. 

Since miRNA functions by regulating its target genes, we may deduce the effects of miRNAs by 
analyzing their regulated networks. To use such a method to elucidate miRNA fiinctions, targets of 
miRNAs should be deduced. Currently, predictions in most target prediction database are based on 
sequence and statistical methods [21]. For example, in TargetScan, seed base pairing, target site 
context, conservation of target site and miRNA, and site accessibility are considered in the prediction 
process [22]. 
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Another method to elucidate miRNA targets is to integrate expression profiles of miRNA and 
mRNA. In the work of Huang et al. [23], a Bayesian-based algorithm, GenMiR++, was developed to 
predict possible targets of 104 miRNAs in humans. They also verified their results by RT-PCR and 
microarray experiments. However, the power of other sequence-based target prediction algorithms was 
not utilized in their work. 

It is also possible to combine sequence-based target prediction and expression-based target 
prediction methods. By integrating expression data into sequence-based predictions, possible false 
positives can be reduced. Previously, miRNA-mRNA interactions were explored with splitting-averaging 
Bayesian networks [24]. In that work, expression profiles of miRNA and mRNA from public 
databases, miRNA target prediction databases, and miRNA sequence information were integrated 
together to discover miRNA-mRNA interaction networks. 

Here, we combined expression profiles of miRNA and mRNA, and three target prediction 
databases, TargetScan, PicTar and miRanda, to obtain confident miRNA-mRNA relationships and 
construct miRNA-regulated protein-protein networks for breast cancer. Furthermore, we explored the 
flinctions of the miRNAs by inspecting the underlying protein interaction networks (PINs) of the 
miRNAs with functional enrichment analysis. This method, as described in Figure 1, was used to 
elucidate the functions of gastric cancer-related miRNAs in our previous work [25]. In that study, a 
gastric cancer-associated miRNA, miR-148a, was identified and validated as being involved in tumor 
proliferation, invasion, migration, and the survival rate of the patients. By using a similar method, we 
aimed to elucidate breast cancer-related miRNA-regulated PINs and their functions. 

Figure 1, Analysis flow chart used in this work. After expression profiles and target 
prediction databases were fetched and preprocessed, they were subjected to the analysis 
process described here and in the "Experimental Section". 
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2. Results and Discussion 

To construct miRNA-regulated PINs, differentially expressed miRNAs and genes from the dataset 
from Farazi et al. [26] were extracted following proper processing of the expression profiles. From our 
selected public miRNA dataset, we found 89 down-regulated miRNAs (93 prior to fold-change 
filtering) and only 1 up-regulated miRNA (Table SI). In gene expression dataset GSE29174, we found 
a total of 1268 down-regulated genes and 587 up-regulated genes before applying the fold change 
filter. There were 726 down-regulated genes (Table S2) and 437 up-regulated genes (Table S3) after 
significantly and differentially expressed genes were filtered by fold change (fold change >2). 

From the results of SAM analysis, we identified some well-known breast cancer-related miRNAs 
(Table SI). For example, miR-214-3p [27] and miR-335-5p [28] have been previously reported to be 
down-regulated in breast cancer. Let-7c was found to be down-regulated in this work, while let-7a, 
another member of the let-7 family, was found to be down-regulated in another work [29]. 
MicroRNAs of the let-7 family were also reportedly down-regulated in several types of cancer [30]. 
We also found that miR-21-5p, the sole up-regulated miRNA in our list, was also previously 
found to be up-regulated [14,31]. However, changes in the expression of most of the miRNAs in our 
down-regulated list have not been reported in the literature. Therefore, we could not rule out the 
possibility that these miRNAs were novel breast cancer-related miRNAs. There are also some 
well-known miRNAs not presented in our list (for such a list, one may see [32-34]). The reason that 
some known miRNAs, for example, miR-19a, miR-155 and miR-205, did not show up in our result 
might be that we used a very stringent threshold (described in Experimental Section) when selecting 
differentially expressed miRNAs for PIN construction. 

Since the miRNAs of the miRNA-regulated PINs were differentially expressed between normal and 
tumor tissues, and we identified some cancer-related functions in our functional enrichment analysis, 
the miRNAs may potentially be useful diagnostic markers for breast cancer. To verify this, we applied 
ROC curve analysis on the miRNA expression profile that was not used in constructing the 
miRNA-regulated PINs. Notably, our results (Figures 2 and SI and Table S4) showed that let-7c 
(Figure 2), miR-497-5p, miR-125b-5p, and some other miRNAs of miRNA-regulated PINs, performed 
well when used as breast cancer diagnostic markers. 

Figure 2. Receiver operating characteristic (ROC) curve of let-7c from our miRNA array 
dataset. For ROC curves of other miRNAs, see Figure S 1 . 
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Following elucidation of differentially expressed miRNAs and genes, miRNA-regulated PINs could 
then be constructed. We identified and constructed partial networks, containing the miRNA and its 
direct target, with the differentially expressed miRNAs and genes as described in the "Experimental 
Section". We then extended the network by appending known interactions from the HPRD 
database. Finally, 18 miRNA-regulated PINs were constructed by the steps described above (Figure 3, 
Figures S2-S13, and Table S5). 

Figure 3. The let-7c-regulated protein interaction network (PIN). This is one of the 18 
miRNA-regulated PINs constructed in this work. Figures of all other miRNA-regulated 
PINs are displayed in Supplementary Figures S2-S13. 
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After construction of the 1 8 PINs was completed, we observed that the sizes of the PINs were not 
similar: some of the miRNAs seemed to regulate larger sized PINs, while other miRNAs affected only 
a small number of genes. Small miRNA-regulated PINs may be caused by the strict g-value threshold 
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set during SAM analysis, tlie processing steps performed on the target prediction databases discussed 
previously, and possibly by lack of protein-protein interaction data for some proteins in the HPRD. 
Although the HPRD may be considered the most comprehensive source of protein-protein interaction 
data [35], some proteins may not have been considered and researched by other investigators, and 
therefore, interaction data for those proteins would not be included in the HPRD. However, it may be 
true that some of the miRNA-regulated PINs were small in breast cancer, since the construction of the 
PINs were based on differentially expressed miRNAs and genes between normal tissues and tumor 
samples, and those miRNAs with small PINs may not be as important as others with larger PINs. 

To elucidate the functions of a miRNA with a regulated PIN, GO enrichment analysis was applied 
to the miRNA-regulated PINs. We did not consider the miRNAs with <5 genes in their regulated PINs, 
and some of the miRNA-regulated PINs had no enriched functions using the defined threshold, 
FDR < 0.0001. To exclude GO terms that describe a broad range of concepts, we only included high 
level GO terms, i.e., larger than 5. 

The results of the GO enrichment analysis for let-7c related to cancer are listed in Table 1. 
(Results of all miRNA-regulated PINs were in Table S6). We defined a GO term as cancer-related if a 
GO term contained "cell proliferation", "cell death", "apoptosis", "signaling", "microtubule", and "actin". 
We noted that 7 miRNAs had enriched GO terms related to apoptosis, cell death, and cell proliferation, 
i.e., miR-520d-3p, miR-497-5p, miR-125b-5p, miR-21-5p, miR-31-5p, let-7c, and miR-125-5p. 
Further, some miRNA-regulated PINs may have functions other than cell survival. For example, the 
nerve growth factor receptor pathway was enriched in miR-regulated PINs of miR-520d-3p, miR-497-5p, 
miR-125a-5p, miR-125b-5p, and miR-31-5p, and the epidermal growth factor receptor pathway was 
enriched in miR-regulated PINs of miR-520d-3p, miR-21-5p, and miR-497-5p. Most of the miRNAs 
had been previously described and were known to be implicated in breast cancer. Let-7c was not only 
likely to be down-regulated in breast cancer [29], but was also found to be a tumor suppressor in 
prostate cancer [16]. Another reported tumor suppressor was miR-125b-5p, which was found to be 
down-regulated in breast tumor tissue [18], and this finding was consistent with our fijnctional 
enrichment results. The miR-125a-5p-regulated PIN was found to be able to inhibit apoptosis and 
regulate epithelial cell proliferation, and has been reported to repress cell growth [36]. 



Table 1. Selected enriched functions of let-7c. Member genes of the let-7c-regulated 
network annotated with corresponding enriched functions are listed. 



MIMAT0000064(hsa-let-7c) 


GO term 


Genes 


Adj. value 


GO:0043067, Regulation of programmed 
cell death 


RRM2B, ACVRl, TP53, RASAl, TGFBRl, PSMA3, BIRC5, 
ACTN2, H0XA13, IRS2, FASTK, VAVl, PSMB6, BCL2, 
CDKl, HDACl, SOXIO, TIAl, AKTl, AURKB 


3.43 X 10-8 


GO:0043069, Negative regulation of 
programmed cell death 


RRM2B, ACVRl, TP53, RASAl, TGFBRl, BIRC5, IRS2, 

BCL2, CDKl, HDACl, SOXIO, AKTl, AURKB 


6.58 X 10"* 


GO:0060548, Negative regulation of cell death 


RRM2B, ACVRl, TP53, RASAl, TGFBRl, B1RC5, IRS2, 
BCL2, CDKl, HDACl, SOXIO, AKTl, AURKB 


8.92 X 10"* 


GO:0ni5630. Microtubule cytoskclcton 


INCENP, SNTB2, SEPTl, TACCl, B1RC5, RACGAPl, P1N4, 
CDCA8, CDKl, I'Hl'l, AKIT, AURKB, NINL, CCDC85B 


5.15 X IQ-' 
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We also found that some breast cancer-specific functions were enriched in our results. For example, 
in the miR-497-5p-reguIated PIN, the term "androgen receptor signaling pathway" was enriched. 
Ahhough it is not clear whether androgens are related to breast cancer, androgen receptors are known 
to be up-regulated in breast cancer and related to node invasiveness [37]. 

To further verify the results of the enriched cancer-related functions, we used GOBO for survival 
analysis. Our hypothesis is that expression of genes annotated with enriched cancer-related terms may 
be related to survival outcome of patients. With exception to some cell death/proliferation-related 
terms, it is already known that some pathways or functions are also related to clinical outcomes. For 
example, cell proliferation-related GO terms have a high probability of affecting survival of cancer 
tissues, and patient outcome may worsen if cancer tissue survives. In addition, some signaling 
pathways were known to enhance invasiveness, migration abilities, or were associated with reduced 
patient survival. For example, the BMP signaling pathway is known to confer various tumor cells with 
enhanced migration and invasion abilities [38], and nerve growth factor receptor (NGFR) was found to 
be associated with overall survival of breast cancer [38]. Furthermore, Toll-like receptor 4 has been 
reported to promote adhesion and invasive migration in breast cancer [39]. Finally, the cytoskeleton 
plays an important role in regulating cell motility in all cells. Actin filaments are known to participate 
in the invasive migration of cancer cells [40]. Since some of these functions were present in our 
enriched terms, we wished to test if the expression of a gene set annotated to the cancer-related 
enriched terms in the PIN would be related to clinical outcome of patients. 

Figure 4. Validation of let-7c result. Gene expression-based outcome for breast cancer 
online (GOBO) survival analysis of let-7c-regulated PIN members marked with the 
following functions: (A) Microtubule cytoskeleton; (B) Negative regulation of 
programmed cell death; and (C) negative regulation of cell death. Red: samples with high 
expression of selected gene set (PIN members); grey: samples with low expression of 
selected gene set (PIN members). 
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As shown in Figure 4, Figure S14 and Table S7, only some of the enriched terms were significantly 
associated with clinical outcome. This may be because the changes in these key genes occurred at the 
protein level, such as in protein expression or even post-translational modification; therefore, mRNA 
expression-centric tools like GOBO cannot explore association of such genes to clinical outcomes. 
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Alternatively, it is possible the miRNA did not regulate the whole pathway, or the miRNA did not 
target the key part of the pathway directly, and thus the clinical outcome of gene sets of the enriched 
GO terms in this condition cannot be determined. However, some functions associated with clinical 
outcomes were observed. For example, proteins annotated with the terms "microtubule cj^oskeleton", 
"negative regulation of programmed cell death", and "negative regulation of cell death" in the 
let-7c-regulated PIN were related to 10 year survival rate of patients, as reported by GOBO (Figure 4). 
Also, the enriched term "regulation of epithelial cell proliferation" for both miR-125a-5p and 
miR-125b-5p were found to be associated with the 10-year survival rate of patients. Therefore, these 
results further supported the GO enrichment analysis discussed previously. 

3. Experimental Section 

3.1. miRNA Microarray Experiments 

We performed a miRNA microarray to obtain the expression profiles for receiver operating 
characteristic (ROC) curve analysis. This dataset was deposited in Gene Expression Omnibus (GEO, 
http://www.ncbi.nlm.nih.gov/geo/, Series Accession: GSE45666). In total, there were 15 normal 
samples and 101 tumor samples in the expression pro tile. Detailed pathophysiogical characteistics of 
these samples were in Table S8 and GSE45666. All human tissue samples collected from breast cancer 
patients were approved and human subject confidentiality was protected by the Institute Review Board 
of National Taiwan University Hospital (IRB, 2007121 IR). 

Total RNA was extracted from tissues collected trom the patients using Trizol*" Reagent (Invitrogen, 
Carlsbad, CA, USA, USA) according to the manufacturer's protocol. Purified RNA was quantified at 
OD260nm by using a ND-1000 spectrophotometer (Nanodrop Technology, Wilmington, DE, USA) 
and qualitated by using a Bioanalyzer 2100 with the RNA 6000 Nano LabChip kit (Agilent Technologies, 
Santa Clara, CA, USA). 

After RNA extraction, 100 ng of total RNA was dephosphorylated and labeled with pCp-Cy3 using 
the Agilent miRNA Complete Labeling and Hyb Kit in conjunction with the microRNA Spike-In kit 
(Agilent Technologies, Santa Clara, CA, USA). Briefly, 2X hybridization buffer (Agilent 
Technologies) was added to the labeled mixture to a final volume of 45 |a.L. The mixture was heated 
for 5 min at 100 °C and immediately cooled to 0 °C. Each 45 \\L sample was hybridized onto an 
Agilent human miRNA Microarray Release 12.0, 8 x 15 K (Agilent Technologies) at 55 °C for 20 h. 
After hybridization, slides were washed for 5 min in Gene Expression Wash Buffer 1 at room 
temperature, then for 5 min in Gene Expression Wash Buffer 2 at 37 °C. Slides were scanned on an 
Agilent microarray scanner (Agilent Technologies, model G2505C) at 100% and 5% sensitivity 
settings. Feature Extraction (Agilent Technologies) software version 10.7.3.1 (Agilent Technologies, 
Santa Clara, CA, USA) was used for image analysis. 

3.2. mRNA Expression Profiles 

For miRNA-regulated protein-protein interaction network construction, the mRNA expression 
profile was fetched from GEO (Series Accession: GSE29174). This dataset was produced by 
Farazi et al. [26], which was the only dataset publicly available with large size of tumor samples and 
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reasonably sized normal samples. In total, 161 clinical samples were collected from breast cancer 
patients by biopsy: 110 invasive ductal carcinoma (IDC), 11 normal, 17 ductal carcinoma in situ 
(DCIS), 1 mucinous A, 8 atypical medullary, 4 apocrine, 8 metaplastic, and 2 adenoid, as classified by 
Farazi et al. Thel 10 IDC samples were classified as the tumor group and the 1 1 normal samples were 
classified as the normal group in this study. 

3.3. miRNA Expression Profiles 

The miRNA expression profile used in this work was fetched from Table S4A of the work of Farazi 
et al. [26]. There were 189 samples in this dataset, with 6 of them from cell lines, and another 183 
samples from patient tissues. In the 183 clinical samples collected from breast cancer patients by 
biopsy in this dataset, there were 128 IDC, 1 1 normal, 18 DCIS, 1 mucinous A, 8 atj^ical medullary, 4 
apocrine, 9 metaplastic, and 2 adenoid cases, as classified by Farazi et al. The 128 IDC samples were 
classified as the tumor group and the 1 1 normal samples as were classified as the normal group. 

3.4. Data Analysis 

The overall workflow design was similar to the previous work of Tseng et al. [25] (see Figure 1). 
However, we applied the workflow on breast cancer expression profiles instead of gastric cancer as in 
the work of Tseng et al. Additionally, we used 3 miRNA target prediction databases, TargetScan 
(v6.0) [22,41,42], PicTar [43,44], and miRanda (release August 2010) [45] here, while only 
TargetScan was used in the previous study. 

To construct the networks, we first elucidated differentially expressed miRNAs and mRNAs 
from published datasets. Both of miRNA (generated with miRNA-Seq technique) and mRNA array 
expression data were produced by Farazi et al. [26] from the same batch of clinical tumor samples. To 
obtain a list of differentially expressed genes and miRNAs between normal and tumor groups, we used 
significance analysis of microarrays (SAM) [46] implemented in R package samr (version 2.0; 
Stanford University, Stanford, CA, USA). We set the false discovery rate (FDR) as <0.0001%, fold 
change as >2.5 for miRNA, and fold change as >1.9 for genes as thresholds to reduce false positives. If 
a gene/miRNA was expressed higher in the normal group compared to the tumor group, we defined 
that gene/miRNA as a down-regulated gene/miRNA, and vice versa. 

Following this, we paired the miRNAs and mRNAs with different expression trends. For example, 
an up-regulated miRNA would be paired with a down-regulated mRNA. If such pairs could be found 
in 2 (or more) of the 3 miRNA target prediction databases, they were then added to their corresponding 
miRNA-regulated network. To further extend the coverage of our network, we incorporated the human 
protein reference database (HPRD) [47], which contains experimentally verified interaction data, into 
our miRNA-regulated PINs. 

Finally, we used gene ontology (GO) enrichment analysis to explore the function of the 
miRNA-regulated PINs. Hypergeomefric tests were used to determine if a GO term was enriched in a 
PIN. In this section, we excluded PINs with less than 5 proteins. We also excluded GO terms with 
levels of less than 5 to avoid non-specific GO terms. Since we tested multiple GO terms on each 
miRNA-regulated PIN, we adjusted the significance of the test with the FDR method developed by 
Benjamini et al. [48]. We used the adjusted value < 0.0001 as our threshold. A GO term would 
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be excluded if its p-wahxQ was larger than 0.0001. Therefore, for each network, we ran the GO 
enrichment analysis, collected the calculated /7-values, and adjusted these values using the methods 
described above. 

3.5. ROC and GOB O Survival Analysis 

After the PINs were constructed, we attempted to verify our results by literature search, ROC, and 
GOBO survival analysis [49]. To determine if the miRNAs we found could serve as classification 
markers for discriminating between normal and tumor samples, we applied ROC analysis on our 
miRNA array data (described in Section 2.1). ROC analysis is usually used to evaluate the efficiency 
of a classifier or a biological marker. R package ROCR [50] (version 1.0.4) was used to plot the ROC 
curve and calculate the area under curve (AUC). The standard error of AUC was then calculated as 
described in the work of Hanlye and McNeil [51]. The /?-value of AUC was thus calculated with 
standard error obtained in the previous step. To further validate if the PINs we found were related to 
cancer, we used survival analysis implemented in GOBO [49] (available at http://co.bmc.lu.se/gobo), 
which provides a large amount of breast cancer gene expression profiles collected from public 
databases with clinical outcome data. In both ROC analysis and GOBO survival analysis, we 
considered our results significant when the /?-value was smaller than 0.05. 

4. Conclusions 

Using integrative analysis of miRNA and mRNA expression profiles, we have identified not only 
breast cancer-related miRNAs and genes, but also putative roles for miRNAs in cancer as elucidated 
fi-om miRNA-regulated PINs constructed in this work. Here, some previously known functions of 
miRNAs were again presented in our results, e.g., the relationship between the miRNAs, let-7c, 
miR-125a-5p, miR-125b-5p, and miR-21-5p, and breast cancer were demonstrated in this research. 
Furthermore, we have identified additional miRNAs and their related fLinctions that have not been 
previously reported or discussed, providing valuable resources for further research in breast cancer. 
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Table SI. Significantly differentially expressed miRNAs found in miRNA dataset in 
Farazi et al. [26]. There are 89 down-regulated miRNAs and 1 up-regulated miRNA in this 
list. 2" values reported by SAM were 0 for all miRNAs in this list. 



nuKBase Accession 


miRNA Name 


Fold Cliange 


miRBase Accession 


miRNA Name 


Fold Change 


MIMAT0004761 


hsa-miR-483-5p 


0.01 


MIMAT0000077 


hsa-miR-22-3p 


0.21 


MIMAT0004552 


hsa-miR-139-3p 


0.01 


MIMAT0000089 


hsa-miR-3 1 -5p 


0.21 


MIMAT0000738 


hsa-miR-383 


0.02 


MIMAT0004612 


hsa-miR- 1 86-3p 


0.21 


MIMAT0002856 


hsa-miR-520d-3p 


0.02 


MIMAT0004592 


hsa-miR- 1 25b- 1 -3p 


0.22 


MIMAT0002811 


hsa-miR-202-3p 


0.03 


MIMAT0001639 


hsa-miR-409-3p 


0.22 


MIMAT0002177 


hsa-miR-486-5p 


0.04 


MIMAT0015032 


hsa-miR-3 1 5 8-3p 


0.22 


MIMAT0022721 


hsa-miR-1247-3p 


0.05 


MIMAT0004496 


hsa-miR-23a-5p 


0.22 


MIMAT0002175 


hsa-miR-485-5p 


0.06 


MIMAT0000690 


hsa-miR-296-5p 


0.22 


MIMAT0000265 


hsa-miR-204-5p 


0.07 


MIMAT0000731 


hsa-miR-3 78a-5p 


0.23 


MIMAT0000752 


hsa-miR-328 


0.07 


MIMAT0000448 


hsa-miR- 1 36-5p 


0.23 


MIMAT0000421 


hsa-miR-122-5p 


0.07 


MIMAT0004796 


hsa-miR-576-3p 


0.23 


MIMAT0000447 


hsa-miR-134 


0.08 


MIMAT0010133 


hsa-miR-2110 


0.23 


MIMAT0000722 


hsa-miR-370 


0.09 


MIMAT0004951 


hsa-miR-887 


0.23 


MIMAT0004513 


hsa-miR-101-5p 


0.09 


MIMAT0003239 


hsa-miR-574-3p 


0.25 


MIMAT0000446 


hsa-miR-127-3p 


0.10 


MIMAT0005901 


hsa-miR- 1249 


0.25 


MIMAT0000097 


hsa-niiR-99a-5p 


0.10 


MIMAT0000510 


hsa-miR-320a 


0.26 


MIMAT0004566 


hsa-miR-218-2-3p 


0.10 


MIMAT0002172 


hsa-miR-376b 


0.26 


MIMAT0000729 


hsa-niiR-376a-3p 


0.11 


MIMAT0000250 


hsa-miR- 139-5p 


0.27 


MIMAT0009197 


hsa-miR-205-3p 


0.11 


MIMAT0005825 


hsa-miR- 1180 


0.27 


MIMAT0004615 


hsa-niiR-195-3p 


0.11 


MIMAT0000437 


hsa-miR- 1 45-5p 


0.28 


MIMAT0005899 


hsa-miR-1247-5p 


0.11 


MIMAT0004601 


hsa-miR- 145-3p 


0.28 


MIMAT0000720 


hsa-miR-376c 


0.12 


MIMAT0003322 


hsa-miR-652-3p 


0.28 


MIMAT0000762 


hsa-miR-324-3p 


0.12 


MIMAT0000756 


hsa-miR-326 


0.28 


MIMAT0004679 


hsa-miR-296-3p 


0.12 


MIMAT0000098 


hsa-miR- 1 00- 5p 


0.29 


MIMAT0004614 


hsa-miR- 1 93a-5p 


0.12 


MIMAT0003296 


hsa-miR-627 


0.29 


MIMAT0003880 


hsa-niiR-67 1 -5p 


0.12 


MIMAT0002820 


hsa-miR-497-5p 


0.31 


MIMAT0004795 


hsa-miR-574-5p 


0.12 


MIMAT0004507 


hsa-miR-92a- 1 -5p 


0.31 


MIMAT0004599 


hsa-miR- 143-5p 


0.13 


MIMAT0000271 


hsa-miR-214-3p 


0.32 


MIMAT0000423 


hsa-miR- 1 25b-5p 


0.13 


MIMAT0004702 


hsa-miR-339-3p 


0.33 


MIMAT0004957 


hsa-miR-760 


0.13 


MIMAT0004611 


hsa-miR- 185-3p 


0.33 


MIMAT0004911 


hsa-miR-874 


0.14 


MIMAT0000064 


hsa-let-7c 


0.34 


MIMAT0004603 


hsa-miR- 1 25b-2-3p 


0.15 


MIMAT0004673 


hsa-miR-29c-5p 


0.35 


MIMAT0004952 


hsa-miR-665 


0.15 


MIMAT0000733 


hsa-miR-379-5p 


0.35 


MIMAT0018205 


hsa-miR-3928 


0.15 


MIMAT0004594 


hsa-miR- 132-5p 


0.35 


MIMAT0004767 


hsa-miR- 193b-5p 


0.15 


MIMAT0000765 


hsa-miR-335-5p 


0.35 


MIMAT0002861 


hsa-miR-518e-3p 


0.15 


MIMAT0002819 


hsa-miR- 193b-3p 


0.36 


MIMAT0004604 


hsa-miR- 127-5p 


0.16 


MIMAT0000088 


hsa-miR-30a-3p 


0.36 


MIMAT0002807 


hsa-miR-491-5p 


0.16 


MIMAT0005951 


hsa-miR- 1307-3p 


0.36 


MIMAT0004689 


hsa-miR-377-5p 


0.16 


MIMAT0004597 


hsa-miR- 140-3p 


0.37 


MIMAT0004762 


hsa-miR-486-3p 


0.16 


MIMAT0004556 


hsa-miR- 10b-3p 


0.37 


MIMAT0000732 


hsa-miR-378a-3p 


0.17 


MIMAT0000272 


hsa-miR-215 


0.37 


MIMAT0017981 


hsa-miR-3605-5p 


0.18 


MIMAT0004511 


hsa-miR-99a-3p 


0.37 


MIMAT0004605 


hsa-miR-129-2-3p 


0.19 


MIMAT0000443 


hsa-miR-125a-5p 


0.38 


MIMAT0006789 


hsa-miR-1468 


0.20 


MIMAT0004482 


hsa-Iet-7b-3p 


0.38 


MIMAT0000737 


hsa-miR-382-5p 


0.21 


MIMAT0000076 


hsa-miR-21-5p 


6.58 
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Table S2. Down-regulated genes found in dataset GSE29174. There are 726 down-regulated 
genes in this list, g- values reported by SAM were 0 for all genes in this list. 



NCBI gene ID 


Gene Symbol 


Fold Change 


NCBI gene ID 


Gene Symbol 


Fold Change 


2949 


GSTM5 


0.06 


619373 


MB0AT4 


0.17 


10894 


LYVEl 


0.06 


130399 


ACVRIC 


0.17 


5950 


RBP4 


0.07 


1646 


AKR1C2 


0.17 


762 


CA4 


0.09 


80763 


C12orf39 


0.17 


54997 


TESC 


0.09 


2159 


FIO 


0.18 


3489 


IGFBP6 


0.09 


84889 


SLC7A3 


0.18 


3952 


LEP 


0.09 


1308 


C0L17A1 


0.18 


213 


ALB 


0.09 


83699 


SH3BGRL2 


0.18 


3131 


HLF 


0.10 


84417 


C2orf40 


0.18 


4023 


LPL 


0.10 


4081 


MAB21L1 


0.18 


10633 


RASLIOA 


0.11 


3484 


IGFBPl 


0.18 


364 


AQP7 


0.11 


5239 


PGM5 


0.19 


1908 


EDN3 


0.11 


4969 


OGN 


0.19 


1811 


SLC26A3 


0.11 


2719 


GPC3 


0.19 


91851 


CHRDLl 


0.11 


116362 


RBP7 


0.19 


729359 


PLES[4 


0.13 


948 


CD36 


0.19 


1149 


CIDEA 


0.13 


5764 


PTN 


0.19 


5959 


RDH5 


0.13 


3043 


HBB 


0.19 


5348 


FXYDl 


0.14 


56920 


SEMA3G 


0.20 


5346 


PLINl 


0.14 


94274 


PPP1R14A 


0.20 


10249 


GLYAT 


0.14 


57447 


NDRG2 


0.20 


158800 


RHOXFl 


0.14 


84795 


PYR0XD2 


0.20 


221476 


PI16 


0.14 


84649 


DGAT2 


0.20 


3040 


HBA2 


0.14 


2690 


GHR 


0.20 


6939 


TCF15 


0.14 


22802 


CLCA4 


0.20 


79645 


EFCABl 


0.14 


5179 


PENK 


0.20 


80343 


SEL1L2 


0.14 


6663 


SOXIO 


0.20 


9413 


FAM189A2 


0.15 


6649 


S0D3 


0.21 


26289 


AK5 


0.15 


54922 


RASIPl 


0.21 


25891 


PAMRl 


0.15 


8406 


SRPX 


0.21 


3679 


ITGA7 


0.15 


1446 


CSNISI 


0.21 


1264 


CNNl 


0.15 


7123 


CLEC3B 


0.22 


92304 


SCGB3A1 


0.15 


9647 


PPMIF 


0.22 


2167 


FABP4 


0.15 


1842 


ECM2 


0.22 


23285 


KIAA1107 


0.15 


3909 


LAMA3 


0.22 


7145 


TNSl 


0.16 


8639 


A0C3 


0.23 


4881 


NPRl 


0.16 


2934 


GSN 


0.23 


1028 


CDKNIC 


0.16 


9370 


ADIPOQ 


0.23 


1036 


CDOl 


0.16 


3202 


HOXA5 


0.23 


130271 


PLEKHH2 


0.16 


9452 


ITM2A 


0.23 


8736 


MYOMl 


0.16 


6290 


SAA3P 


0.23 


8908 


GYG2 


0.16 


4604 


MYBPCl 


0.23 
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Table S2, Cont. 



NCBI gene ID 


Gene Symbol 


Fold Change 


NCBI gene ID 


Gene Symbol 


Fold Change 


79785 


RERGL 


0.16 


1128 


CHRMl 


0.23 


221091 


LRRN4CL 


0.17 


83878 


USHBPl 


0.24 


3991 


LIFE 


0.17 


63970 


TP53AIP1 


0.24 


27175 


TUBG2 


0.24 


79192 


IRXl 


0.28 


1346 


C0X7A1 


0.24 


3400 


ID4 


0.28 


6376 


CX3CL1 


0.24 


57519 


STARD9 


0.29 


50486 


G0S2 


0.24 


57666 


FBRSLl 


0.29 


6285 


SI GOB 


0.24 


3590 


ILllRA 


0.29 


443 


ASPA 


0.24 


57664 


PLEKHA4 


0.29 


947 


CD34 


0.25 


197257 


LDHD 


0.29 


84632 


AFAP1L2 


0.25 


66036 


MTMR9 


0.29 


3866 


KRT15 


0.25 


2321 


FLTl 


0.29 


147463 


ANKRD29 


0.25 


126 


ADHIC 


0.29 


2878 


GPX3 


0.25 


1363 


CPE 


0.29 


7079 


TEV[P4 


0.25 


56131 


X^ I X XX^ A 

PCDHB4 


0.29 


54345 


SOX18 


0.25 


22915 


MMRNl 


0.29 


51277 


DNAJC27 


0.25 


7069 


rnx XXX Xfc 

THRSP 


0.29 


84870 


RSP03 


0.25 


57161 


PELI2 


0.30 


55323 


LARP6 


0.25 


770 


CAll 


0.30 


6387 


CXCL12 


0.25 


53342 


IL17D 


0.30 


137835 


TMEM71 


0.25 


79987 


SVEPl 


0.30 


5212 


VIT 


0.25 


857 


CAVl 


0.30 


26577 


PC0LCE2 


0.25 


222166 


C7orf41 


0.30 


845 


CASQ2 


0.25 


27190 


IL17B 


0.30 


6422 


SFRPl 


0.25 


116159 


CYYRl 


0.30 


10351 


ABCA8 


0.26 


4487 


MSXl 


0.30 


10840 


ALDHILI 


0.26 


9068 


ANGPTLl 


0.30 


65983 


GRAMD3 


0.26 


10411 


RAPGEF3 


0.30 


84327 


ZBED3 


0.26 


3199 


H0XA2 


0.30 


57124 


CD248 


0.26 


2944 


GSTMl 


0.30 


3235 


HOXD9 


0.26 


2920 


CXCL2 


0.30 


2192 


FBLNl 


0.26 


201134 


CEP112 


0.31 


91653 


BOG 


0.26 


220001 


VWCE 


0.31 


4147 


MATN2 


0.26 


83888 


xn 1 1 1 f V x^ ^ 

FGFBP2 


0.31 


126669 


SHE 


0.27 


6366 


CCL21 


0.31 


2788 


GNG7 


0.27 


6711 


SPTBNl 


0.31 


129804 


FBLN7 


0.27 


85378 


TUBGCP6 


0.31 


270 


AMPDl 


0.27 


26040 


SETBPl 


0.31 


79656 


BENDS 


0.27 


4692 


NDN 


0.31 


58503 


PROLl 


0.27 


25890 


ABI3BP 


0.31 


3316 


HSPB2 


0.27 


23531 


MMD 


0.31 


729440 


CCDC61 


0.27 


30846 


EHD2 


0.31 


54438 


GFODl 


0.27 


6196 


RPS6KA2 


0.31 


5243 


ABCBl 


0.27 


2009 


EMLl 


0.31 
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Table S2, Cont. 



NCBI gene ID 


Gene Symbol 


Fold Change 


NCBI gene ID 


Gene Symbol 


Fold Change 


810 


CALML3 


0.27 


6289 


SAA2 


0.31 


6898 


TAT 


0.27 


345275 


HSD17B13 


0.31 


5648 


MASPl 


0.28 


2701 


GJA4 


0.32 


25999 


CLIPS 


0.28 


112609 


MRAP2 


0.32 


125875 


CLDND2 


0.28 


727 


C5 


0.S2 


7102 


TSPAN7 


0.28 


477 


ATP1A2 


0.32 


1879 


EBFl 


0.28 


9627 


SNCAIP 


0.S2 


23252 


OTUD3 


0.28 


4435 


CITED 1 


0.32 


5493 


PPL 


0.28 


10974 


C10orfll6 


0.32 


83987 


CCDC8 


0.28 


11005 


SPINK5 


0.32 


9073 


CLDN8 


0.28 


80325 


ABTBl 


O.SS 


221981 


THSD7A 


0.28 


221395 


GPR116 


O.SS 


64102 


TNMD 


0.28 


10014 


HDAC5 


O.SS 


137872 


ADHFEl 


O.SS 


1489 


CTFl 


0.37 


27151 


CP AMDS 


O.SS 


35 


ACADS 


0.37 


387923 


SERP2 


O.SS 


3749 


KCNC4 


0.37 


145581 


LRFN5 


O.SS 


140738 


TMEM37 


0.S7 


6263 


RYR3 


O.SS 


2791 


GNGll 


0.S7 


2354 


FOSB 


O.SS 


23604 


DAPK2 


0.S7 


51302 


CYP39A1 


O.SS 


10217 


CTDSPL 


0.37 


4128 


MAOA 


0.S4 


23550 


PSD4 


0.37 


117248 


GALNTL2 


0.S4 


4306 


NRSC2 


0.37 


10268 


RAMPS 


0.S4 


119587 


CPXM2 


0.S7 


7730 


ZNF177 


0.S4 


7942 


TFEB 


0.S7 


10873 


MES 


0.S4 


3815 


KIT 


0.S7 


7461 


CLIP2 


0.S4 


1805 


DPT 


0.37 


7049 


TGFBRS 


0.S4 


23242 


COBL 


0.37 


79901 


CYBRDl 


0.S4 


4313 


MMP2 


0.37 


5152 


PDE9A 


0.S4 


4139 


MARKl 


0.S7 


50805 


IRX4 


0.S4 


9104 


RGN 


0.S7 


8644 


AKRICS 


0.S4 


2329 


FM04 


0.S7 


5915 


RARB 


0.S4 


25802 


LMODl 


0.38 


2770 


GNAIl 


0.S4 


4239 


MFAP4 


0.38 


54996 


2-Mar 


0.35 


10392 


NODI 


0.38 


79791 


FBXOS 1 


0.35 


6794 


STKll 


0.38 


54776 


PPP1R12C 


0.S5 


85458 


DIXDCl 


0.38 


9079 


LDB2 


0.S5 


4123 


MAN2C1 


0.38 


57104 


PNPLA2 


0.S5 


54476 


RNF216 


0.38 


30008 


EFEMP2 


0.S5 


9920 


KBTBDll 


0.38 


91461 


PKDCC 


0.35 


6329 


SCN4A 


0.38 


23368 


PPPIRISB 


0.S5 


10253 


SPRY2 


0.38 


23461 


ABCA5 


0.S5 


1910 


EDNRB 


0.38 


9572 


NRIDI 


0.35 


9249 


DHRS3 


0.38 


23338 


PHF15 


0.35 


22869 


ZNF510 


0.38 
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Table S2, Cont. 



NCBI gene ID 


Gene Symbol 


Fold Change 


NCBI gene ID 


Gene Symbol 


Fold Change 


114800 


CCDC85A 


0.35 


3384 


ICAM2 


0.38 


2550 


GABBRl 


0.35 


8613 


PPAP2B 


0.38 


4638 


MYLK 


0.35 


1950 


EGF 


0.38 


2327 


FM02 


0.35 


55273 


TMEMIOO 


0.38 


139411 


PTCHDl 


0.35 


6297 


SALL2 


0.38 


10391 


y^ x^ y^ ^ x^ 

C0R02B 


0.35 


9365 


KL 


0.38 


25854 


FAM149A 


0.35 


8863 


PER3 


0.38 


55701 


ARHGEF40 


0.36 


8404 


SPARCLl 


0.38 


1759 


DNMl 


0.36 


2202 


EFEMPl 


0.38 


22849 


CPEB3 


0.36 


8369 


HIST1H4G 


0.38 


57716 


PRX 


0.36 


5187 


PERI 


0.39 


1628 


DBF 


0.36 


30815 


<~\ rxn A x ta T A ^ 

ST6GALNAC6 


0.39 


80031 


SEMA6D 


0.36 


256364 


EML3 


0.39 


259217 


HSPA12A 


0.36 


57381 


RHOJ 


0.39 


6909 


TBX2 


0.36 


761 


CA3 


0.39 


1511 


CTSG 


0.36 


83989 


FAM172A 


0.39 


79971 


WLS 


0.36 


1408 


CRY2 


0.39 


90865 


IL33 


0.36 


2281 


FKBPIB 


0.39 


11343 


MGLL 


0.36 


51222 


ZNF219 


0.39 


55800 


SCN3B 


0.36 


54540 


FAM193B 


0.39 


1949 


EFNB3 


0.36 


4053 


LTBP2 


0.39 


284217 


LAMAl 


0.36 


55184 


DZANKl 


0.39 


22927 


HABP4 


0.37 


5740 


PTGIS 


0.39 


23645 


PPP1R15A 


0.39 


84814 


PPAPDC3 


0.42 


342574 


KRT27 


0.39 


79365 


BHLHE41 


0.42 


83543 


AIFIL 


0.39 


316 


AOXl 


0.42 


624 


BDKRB2 


0.39 


23380 


SRGAP2 


0.42 


347 


APOD 


0.39 


84033 


OBSCN 


0.42 


84935 


C13orf33 


0.39 


90353 


CTUl 


0.42 


858 


CAV2 


0.39 


9013 


TAFIC 


0.42 


5138 


PDE2A 


0.40 


474344 


GIMAP6 


0.42 


114928 


GPRASP2 


0.40 


84883 


AIFM2 


0.42 


58190 


CTDSPl 


0.40 


58480 


RHOU 


0.42 


513 


ATP5D 


0.40 


65982 


ZSCAN18 


0.42 


57684 


ZBTB26 


0.40 


666 


BOK 


0.42 


7041 


TGFBIII 


0.40 


79762 


Clorfll5 


0.42 


5787 


PTPRB 


0.40 


525 


ATP6V1B1 


0.42 


7294 


TXK 


0.40 


4675 


NAP1L3 


0.42 


56301 


SLC7A10 


0.40 


3257 


HPSl 


0.43 


55937 


APOM 


0.40 


55781 


RIOK2 


0.43 


6368 


CCL23 


0.40 


63947 


DMRTCl 


0.43 


55020 


TTC38 


0.40 


1969 


EPHA2 


0.43 


134265 


AFAPILI 


0.40 


25927 


CNRIPl 


0.43 


4485 


MSTl 


0.40 


57685 


CACHDl 


0.43 
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Table S2, Cont. 



NCBI gene ID 


Gene Symbol 


Fold Change 


NCBI gene ID 


Gene Symbol 


Fold Change 


51559 


NT5DC3 


0.40 


29997 


GLTSCR2 


0.43 


7169 


TPM2 


0.40 


26051 


PPP1R16B 


0.43 


51705 


EMCN 


0.40 


83604 


TMEM47 


0.43 


8938 


BAIAP3 


0.40 


2308 


FOXOl 


0.43 


10365 


KLF2 


0.40 


55225 


RAVER2 


0.43 


59 


ACTA2 


0.40 


54839 


LRRC49 


0.43 


80309 


SPHKAP 


0.40 


122953 


JDP2 


0.43 


3779 


KCNMBl 


0.41 


29775 


CARDIO 


0.43 


10826 


C5orf4 


0.41 


166 


AES 


0.43 


219654 


ZCCHC24 


0.41 


25924 


MYRIP 


0.43 


92162 


TMEM88 


0.41 


2852 


GPER 


0.43 


7450 


VWF 


0.41 


51421 


AM0TL2 


0.43 


10266 


RAMP2 


0.41 


124936 


CYB5D2 


0.43 


25875 


LETMDl 


0.41 


1294 


COL7A1 


0.43 


1938 


EEF2 


0.41 


127435 


PGDN 


0.43 


121551 


BTBDll 


0.41 


84952 


CGNLl 


0.43 


2119 


ETV5 


0.41 


83483 


PLVAP 


0.43 


9696 


CROCC 


0.41 


1958 


EGRl 


0.43 


1031 


CDKN2C 


0.41 


230 


ALDOC 


0.43 


9037 


SEMA5A 


0.41 


65987 


KCTD14 


0.43 


3397 


IDl 


0.41 


4804 


NGFR 


0.44 


84707 


BEX2 


0.41 


64852 


TUTl 


0.44 


57616 


TSHZ3 


0.41 


84253 


GARNL3 


0.44 


1471 


CST3 


0.41 


5866 


RAB31L1 


0.44 


55214 


LEPRELl 


0.41 


10608 


MXD4 


0.44 


3914 


LAMB3 


0.41 


4211 


MEISl 


0.44 


57478 


USP31 


0.41 


83547 


RILP 


0.44 


3783 


KCNN4 


0.41 


9172 


MYOM2 


0.44 


8839 


WISP2 


0.41 


57192 


MCOLNl 


0.44 


1583 


CYPllAl 


0.42 


255877 


BCL6B 


0.44 


10124 


ARL4A 


0.42 


56904 


SH3GLB2 


0.44 


738 


Cllorf2 


0.42 


51285 


RASL12 


0.44 


29800 


ZDHHCl 


0.42 


3425 


IDUA 


0.44 


23135 


KDM6B 


0.44 


402117 


VWC2L 


0.46 


171024 


SYNP02 


0.44 


81490 


PTDSS2 


0.46 


10350 


ABCA9 


0.44 


^ A mm A ^ 

283748 


PLA2G4D 


0.46 


3691 


ITGB4 


0.44 


23523 


CABIN 1 


0.46 


2348 


FOLRl 


0.44 


6146 


RPL22 


0.46 


11145 


PLA2G16 


0.44 


85360 


SYDEl 


0.46 


554 


AVPR2 


0.45 


60468 


BACH2 


0.46 


64072 


CDH23 


0.45 


57451 


0DZ2 


0.46 


80177 


MYCTl 


0.45 


4013 


VWA5A 


0.46 


5957 


RCVRN 


0.45 


339768 


ESPNL 


0.46 


408 


ARRBl 


0.45 


3860 


KRT13 


0.46 



Int. J. Mol. Sci. 2013, 14 



11576 



Table S2, Cont. 



NCBI gene ID 


Gene Symbol 


Fold Change 


NCBI gene ID 


Gene Symbol 


Fold Change 


144699 


FBXL14 


0.45 


7094 


TLNl 


0.46 


83719 


YPEL3 


0.45 


4232 


MEST 


0.46 


22841 


RAB11FIP2 


0.45 


1410 


CRYAB 


0.46 


283927 


NUDT7 


0.45 


57452 


GALNTLl 


0.47 


293 


SLC25A6 


0.45 


63935 


PCIFl 


0.47 


90507 


SCRN2 


0.45 


25873 


RPL36 


0.47 


37 


ACADVL 


0.45 


9812 


KI A AO 141 


0.47 


112744 


IL17F 


0.45 


51665 


ASBl 


0.47 


6709 


SPTANl 


0.45 


64123 


ELTDl 


0.47 


8086 


AAAS 


0.45 


6122 


RPL3 


0.47 


7423 


VEGFB 


0.45 


222962 


SLC29A4 


0.47 


64221 


R0B03 


0.45 


23102 


TBC1D2B 


0.47 


7273 


TTN 


0.45 


3476 


IGBPl 


0.47 


2657 


GDFl 


0.45 


93408 


MYLIO 


0.47 


59271 


C21orf63 


0.45 


5310 


PKDl 


0.47 


132160 


PPMIM 


0.45 


4628 


MYHIO 


0.47 


27244 


SESNl 


0.45 


221935 


SDKl 


0.47 


51310 


SLC22A17 


0.45 


23328 


SASHl 


0.47 


4828 


NMB 


0.45 


8522 


GAS 7 


0.47 


54360 


CYTLl 


0.45 


10023 


FRATl 


0.47 


203245 


NAIFl 


0.45 


7301 


TYR03 


0.47 


23166 


STABl 


0.45 


2767 


GNAll 


0.47 


2121 


EVC 


0.45 


9457 


FHL5 


0.47 


116496 


FAM129A 


0.45 


4094 


MAF 


0.47 


23239 


PHLPPl 


0.45 


65268 


WNK2 


0.47 


51673 


rill 

TPPP3 


0.45 


54585 


T" r — 7 r 1 1 1 1 T i 

LZTFLl 


0.47 


64094 


SMOC2 


0.45 


375449 


MAST4 


0.47 


6383 


SDC2 


0.45 


138311 


FAM69B 


0.47 


2180 


ACSLl 


0.45 


160622 


GRASP 


0.47 


23770 


FKBP8 


0.45 


22837 


COBLLl 


0.47 


55901 


THSDl 


0.46 


51435 


SCARA3 


0.47 


25895 


METTL21B 


0.46 


217 


ALDH2 


0.47 


23731 


C9orf5 


0.46 


6236 


RRAD 


0.47 


126393 


HSPB6 


0.46 


8322 


FZD4 


0.47 


4056 


LTC4S 


0.46 


653275 


CFCIB 


0.47 


79825 


CCDC48 


0.46 


10908 


PNPLA6 


0.47 


10810 


WASF3 


0.46 


57526 


PCDH19 


0.47 


29911 


HOOK2 


0.46 


8424 


BBOXl 


0.47 


583 


BBS2 


0.46 


9905 


SGSM2 


0.48 


28984 


C13orfl5 


0.46 


10435 


CDC42EP2 


0.48 


1465 


CSRPl 


0.46 


23087 


TRIM35 


0.48 


55258 


THNSL2 


0.46 


60314 


C12orfl0 


0.48 


161198 


CLEC14A 


0.46 


1073 


CFL2 


0.48 


3699 


ITffl3 


0.48 


5256 


PHKA2 


0.49 
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Table S2, Cont. 



NCBI gene ID 


Gene Symbol 


Fold Change 


NCBI gene ID 


Gene Symbol 


Fold Change 


92922 


CCDC102A 


0.48 


6237 


RRAS 


0.49 


65057 


ACD 


0.48 


5288 


PIK3C2G 


0.49 


9095 


TBX19 


0.48 


10252 


SPRYl 


0.49 


6441 


£^ rill 

SFTPD 


0.48 


79026 


AHNAK 


0.49 


22846 


VASHl 


0.48 


9693 


RAPGEF2 


0.49 


51066 


C3orf32 


0.48 


51226 


C0PZ2 


0.49 


23179 


RGLl 


0.48 


158326 


FREMl 


0.49 


4664 


NABl 


0.48 


1956 


EGFR 


0.49 


50511 


SYCP3 


0.48 


5360 


PLTP 


0.49 


6430 


SRSF5 


0.48 


290 


ANPEP 


0.49 


11078 


TRIOBP 


0.48 


1756 


DMD 


0.49 


78991 


PCYOXIL 


0.48 


5118 


PCOLCE 


0.49 


6623 


SNCG 


0.48 


56654 


NPDCl 


0.49 


23384 


SPECCIL 


0.48 


9254 


CACNA2D2 


0.49 


53826 


FXYD6 


0.48 


55536 


CDCA7L 


0.49 


9397 


NMT2 


0.48 


124975 


GGT6 


0.49 


6041 


RNASEL 


0.48 


1906 


EDNl 


0.49 


113510 


HELQ 


0.48 


81029 


WNT5B 


0.49 


64788 


LMFl 


0.48 


2646 


GCKR 


0.49 


2217 


FCGRT 


0.48 


9811 


CTIF 


0.50 


79720 


VPS37B 


0.48 


145376 


PPP1R36 


0.50 


6764 


ST5 


0.48 


222865 


TMEM130 


0.50 


252969 


NEIL2 


0.48 


rf~V A rf~v rf~V rf~V 

92999 


ZBTB47 


0.50 


8987 


STBDl 


0.48 


168002 


DACT2 


0.50 


41 


ACCN2 


0.48 


6829 


SUPT5H 


0.50 


7905 


REEP5 


0.48 


9992 


X T" y^"V Txn ^ 

KCNE2 


0.50 


5919 


RARRES2 


0.48 


58509 


C19orf29 


0.50 


10544 


PROCR 


0.48 


79706 


X^X^ X7"X^ XX^ 1 

PRKRIPl 


0.50 


6876 


TAGLN 


0.48 


1153 


CIRBP 


0.50 


8436 


SDPR 


0.49 


9639 


ARHGEFIO 


0.50 


23500 


DAAM2 


0.49 


4054 


LTBP3 


0.50 


130132 


RFTN2 


0.49 


1120 


CHKB 


0.50 


80310 


PDGFD 


0.49 


286046 


XKR6 


0.50 


4215 


MAP3K3 


0.49 


9590 


AKAP12 


0.50 


282775 


OR5J2 


0.49 


64115 


C10orf54 


0.50 


51161 


C3orfl8 


0.49 


2067 


ERCCl 


0.50 


29098 


RANGRF 


0.49 


7507 


XPA 


0.50 


53336 


CPXCRl 


0.49 




CEP 164 


0.50 


9081 


PRY 


0.49 


652 


BMP4 


0.50 


9459 


ARHGEF6 


0.49 


55702 


CCDC94 


0.50 


2995 


GYPC 


0.49 


57613 


KIAA1467 


0.50 


23057 


NMNAT2 


0.49 


28514 


DLLl 


0.50 


4669 


NAGLU 


0.49 


169270 


ZNF596 


0.50 


6452 


SH3BP2 


0.49 


83982 


IFI27L2 


0.50 
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Table S2, Cont. 



JNCBl gene ID 


Gene Symbol 


*ola Change 


JNCBl gene ID 


y~i o ■ „ 1 

Gene Symbol 


Fold Change 


51458 


RHCG 


0.49 


2247 


FGF2 


0.50 


1112 


FOXN3 


0.49 


26248 


OR2K2 


0.50 


29954 


POMT2 


0.49 


84303 


CHCHD6 


0.50 


9612 


NCOR2 


0.49 


3615 


IMPDH2 


0.50 


3198 


HOXAl 


0.49 


1813 


DRD2 


0.50 


5311 


PKD2 


0.49 


80148 


PQLCl 


0.50 


2946 


GSTM2 


0.49 


390081 


OR52E4 


0.50 


2109 


ETFB 


0.49 


352954 


GATS 


0.50 


56062 


KLHL4 


0.49 


90871 


C9orfl23 


0.50 


6915 


TBXA2R 


0.50 


50945 


TBX22 


0.52 


64288 


ZNF323 


0.50 


5204 


PFDN5 


0.52 


5195 


PEX14 


0.50 


5338 


PLD2 


0.52 


84557 


MAP1LC3A 


0.50 


94 


ACVRLl 


0.52 


6164 


RPL34 


0.50 


54039 


PCBP3 


0.52 


8835 


SOCS2 


0.50 


7691 


r^"v 111 -t ^ ^ 

ZNF132 


0.52 


2735 


GLIl 


0.50 


338 


APOB 


0.52 


26022 


TMEM98 


0.50 


84658 


EMR3 


0.52 


3908 


LAMA2 


0.50 


283232 


TMEM80 


0.52 


1825 


DSC3 


0.50 


5430 


P0LR2A 


0.52 


5730 


PTGDS 


0.50 


54623 


PAFl 


0.52 


162515 


SLC16A11 


0.51 


11070 


TMEM115 


0.52 


274 


BESFl 


0.51 


10395 


DLCl 


0.52 


79654 


HECTD3 


0.51 


57140 


RNPEPLl 


0.52 


22863 


ATG14 


0.51 


79781 


IQCAl 


0.52 


25949 


SYF2 


0.51 


1838 


DTNB 


0.52 


84872 


ZC3H10 


0.51 


51386 


EIF3L 


0.52 


23187 


PHLDBl 


0.51 


56919 


DHX33 


0.52 


5434 


POLR2E 


0.51 


57542 


KLHDC5 


0.52 


6181 


RPLP2 


0.51 


3628 


INPPl 


0.52 


6141 


RPL18 


0.51 


4520 


MTFl 


0.52 


84747 


LrNC119B 


0.51 


8547 


FCN3 


0.52 


23399 


CTDNEPl 


0.51 


60401 


EDA2R 


0.52 


599 


BCL2L2 


0.51 


8082 


SSPN 


0.52 


197258 


FUK 


0.51 


80755 


AARSDl 


0.52 


5207 


PFKFBl 


0.51 


710 


SERPINGI 


0.52 


8131 


NPRL3 


0.51 


56246 


MRAP 


0.52 


25839 


C0G4 


0.51 


10555 


A y^ A rxn/^ 

AGPAT2 


0.52 


10816 


SPINT3 


0.51 


949 


SCARBl 


0.52 


60485 


SAVl 


0.51 


23743 


BHMT2 


0.52 


5681 


PSKHl 


0.51 


3910 


LAMA4 


0.52 


80318 


GKAPl 


0.51 


60370 


AVPll 


0.52 


57088 


PLSCR4 


0.51 


5021 


OXTR 


0.52 


93129 


0RAI3 


0.51 


55997 


CFCl 


0.52 


5829 


PXN 


0.51 


23144 


ZC3H3 


0.52 
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Table S2, Cont. 



JNCBl gene ID 


Gene Symbol 


*ola Change 


JNCBl gene ID 


Gene Symbol 


"17' „ 1 J 

*old Change 


56776 


FMN2 


0.51 


150709 


ANKAR 


0.52 


85456 


TNKSIBPI 


0.51 


6591 


SNAI2 


0.52 


283 


ANG 


0.51 


10129 


FRY 


0.52 


7035 


TFPI 


0.51 


5166 


PDK4 


0.52 


51232 


CRIMl 


0.51 


146433 


1L34 


0.52 


112616 


CMTM7 


0.51 


118812 


M0RN4 


0.53 


22981 


NINL 


0.51 


10516 


FBLN5 


0.53 


8727 


CTNNALl 


0.51 


9463 


PICKl 


0.53 


9902 


MRC2 


0.51 


127495 


LRRC39 


0.53 


10900 


RUNDC3A 


0.51 


7753 


ZNF202 


0.53 


51299 


NRNl 


0.51 


79827 


CLMP 


0.53 


79632 


FAM184A 


0.52 


203260 


CCDC107 


0.53 


80820 


EEPDl 


0.52 


83657 


DYNLRB2 


0.53 


Table SI. Up-regulated genes found in dataset GSE29174. There are 437 


up-regulated 


genes in this list, values reported by SAM were 0 for all genes in this list. 




NCBI gene ID 


Gene Symbol 


Fold Change 


NCBI gene ID 


Gene Symbol 


Fold Change 


1300 


COLlOAl 


42.74 


54443 


ANLN 


5.79 


3007 


HISTIHID 


29.72 


6710 


SPTB 


5.71 


8366 


HIST1H4B 


25.58 


7272 


TTK 


5.64 


6286 


SI OOP 


25.19 


10635 


RAD51AP1 


5.49 


1301 


COLllAl 


24.72 


4069 


LYZ 


5.37 


3627 


CXCLIO 


17.83 


55183 


RIFl 


5.34 


4283 


CXCL9 


15.88 


891 


CCNBl 


5.34 


1387 


CREBBP 


12.83 


91543 


RSAD2 


5.31 


27299 


ADAMDECl 


12.78 


81610 


FAM83D 


5.24 


54986 


ULK4 


12.46 


64581 


CLEC7A 


5.10 


55771 


PRRll 


12.02 


10051 


SMC4 


5.02 


54790 


TET2 


11.25 


4085 


MAD2L1 


4.96 


6241 


RRM2 


10.60 


55872 


PBK 


4.83 


3433 


IFIT2 


10.49 


991 


CDC20 


4.82 


6999 


TD02 


9.73 


9221 


NOLCl 


4.74 


1656 


DDX6 


9.72 


2124 


EV12B 


4.66 


55088 


C10orfll8 


9.37 


375248 


ANKRD36 


4.66 


9648 


GCC2 


9.24 


1164 


CKS2 


4.64 


6696 


SPPl 


8.92 


1230 


CCRl 


4.62 


2803 


G0LGA4 


8.57 


890 


CCNA2 


4.56 


83540 


NUF2 


7.73 


127933 


UHMKl 


4.49 


10112 


KIF20A 


7.66 


10274 


STAGl 


4.45 


9833 


MELK 


7.59 


597 


BCL2A1 


4.43 


55165 


CEP55 


7.50 


55355 


HJURP 


4.41 


10142 


AKAP9 


7.44 


54210 


TREMl 


4.36 


9447 


AIM2 


7.42 


253558 


LCLATl 


4.26 
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Table S3. Cont. 



NCBl gene ID 


Gene Symbol 


Fold Change 


JNCBl gene ID 


Gene Symbol 


Fold Change 


2706 


GJB2 


7.33 


1033 


CDKN3 


4.24 


6498 


SKK. 


7.13 


79801 


SHCBPl 


4.23 


219285 


SAMD9L 


7.06 


126731 


Clorf96 


4.21 


10261 


IGSF6 


7.01 


6772 


STATl 


4.20 


2335 


FNl 


6.95 


55729 


ATF7IP 


4.14 


699 


BUBl 


6.75 


6713 


SQLE 


4.14 


1058 


CENPA 


6.75 


157570 


ESC02 


4.10 


332 


BIRC5 


6.73 


79871 


RPAP2 


4.09 


51203 


NUSAPl 


6.59 


9493 


KIF23 


4.09 


259266 


ASPM 


6.54 


4751 


NEK2 


4.05 


1063 


CENPF 


6.49 


10631 


POSTN 


4.03 


165918 


RNF168 


6.44 


23515 


M0RC3 


4.02 


9232 


PTTGl 


6.34 


7153 


T0P2A 


4.02 


5996 


RGSl 


6.07 


10403 


NDC80 


4.00 


29089 


UBE2T 


5.96 


10915 


TCERGl 


3.99 


22974 


TPX2 


5.94 


57650 


KIAA1524 


3.99 


4321 


MMP12 


5.91 


23049 


SMGl 


3.93 


983 


CDKl 


5.89 


80231 


CXorf21 


3.87 


85444 


LRRCCl 


5.87 


5111 


PCNA 


3.86 


29121 


CLEC2D 


3.83 


79682 


MLFIIP 


3.11 


4090 


SMAD5 


3.80 


29123 


ANKRDll 


3.09 


2123 


EVI2A 


3.80 


5429 


POLH 


3.09 


57695 


USP37 


3.79 


701 


BUB IB 


3.07 


133418 


EMB 


3.76 


200030 


NBPFU 


3.06 


4131 


MAP IB 


3.76 


55677 


IWSl 


3.06 


9787 


DLGAP5 


3.75 


160418 


TMTC3 


3.04 


9768 


KIAAOlOl 


3.74 


9147 


NEMF 


3.04 


54625 


PARP14 


3.73 


11320 


MGAT4A 


3.04 


2215 


FCGR3B 


3.71 


5238 


PGM3 


3.03 


9134 


CCNE2 


3.70 


2820 


GPD2 


3.02 


3117 


HLA-DQAl 


3.68 


388886 


FAM211B 


3.01 


10380 


BPNTl 


3.67 


7852 


CXCR4 


3.00 


79056 


PRRG4 


3.63 


57082 


CASC5 


2.99 


10673 


TNFSF13B 


3.63 


22926 


ATF6 


2.98 


8467 


SMARCA5 


3.61 


7594 


ZNF43 


2.98 


115908 


CTHRCl 


3.61 


968 


CD68 


2.97 


3428 


IFI16 


3.61 


7171 


TPM4 


2.96 


1520 


CTSS 


3.61 


11004 


KIF2C 


2.96 


10797 


MTHFD2 


3.57 


10808 


HSPHl 


2.95 


55681 


SCYL2 


3.57 


84909 


C9orf3 


2.94 


9749 


PHACTR2 


3.57 


1894 


ECT2 


2.93 


94240 


EPSTIl 


3.56 


1629 


DBT 


2.92 


64151 


NCAPG 


3.51 


116969 


ART5 


2.90 


25879 


DCAF13 


3.51 


3227 


HOXCll 


2.88 
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Table S3. Cont. 



NCBI gene ID 


Gene Symbol 


Fold Change 


NCBI gene ID 


Gene Symbol 


Fold Change 


116064 


LRRC58 


3.47 


3149 


HMGB3 


2.87 


29899 


GPSM2 


3.47 


10437 


IFI30 


2.87 


135114 


HINTS 


3.45 


57489 


0DF2L 


2.87 


27333 


G0LIM4 


3.43 


2151 


F2RL2 


2.86 


55839 


CENPN 


3.43 


23215 


PRRC2C 


2.85 


23213 


SULFl 


3.41 


128710 


C20orf94 


2.85 


81671 


VMPl 


3.39 


23594 


0RC6 


2.84 


9889 


ZBED4 


3.36 


5205 


ATP8B1 


2.83 


3092 


HIPl 


3.34 


51430 


Clorf9 


2.80 


51512 


GTSEl 


3.34 


57405 


SPC25 


2.80 


92797 


HELB 


3.34 


112401 


BIRC8 


2.80 


51426 


POLK 


3.30 


3606 


IL18 


2.80 


5611 


DNAJC3 


3.30 


115362 


GBP5 


2.80 


6596 


HLTF 


3.28 


50515 


CHSTl 1 


2.79 


9910 


RABGAPIL 


3.25 


83461 


CDCA3 


2.79 


528 


ATP6V1C1 


3.23 


10744 


PTTG2 


2.78 


3833 


KIFCl 


3.23 


51765 


MST4 


2.77 


197131 


UBRl 


3.20 


10926 


DBF4 


2.76 


29923 


HILPDA 


3.20 


27125 


AFF4 


2.75 


28998 


MRPL13 


3.19 


10615 


SPAG5 


2.75 


58527 


C6orfll5 


3.19 


55143 


CDCA8 


2.74 


79000 


Clorfl35 


3.19 


51602 


NOP58 


2.74 


9857 


CEP350 


3.18 


51478 


HSD17B7 


2.73 


84296 


GINS4 


3.18 


2209 


FCGRIA 


2.73 


81034 


SLC25A32 


3.15 


9958 


USP15 


2.72 


55723 


ASF IB 


3.14 


5469 


MEDl 


2.72 


7110 


TMFl 


3.14 


8813 


DPMI 


2.70 


84081 


NSRPl 


3.14 


6731 


SRP72 


2.70 


23075 


SWAP70 


3.12 


9991 


PTBP3 


2.70 


6726 


SRP9 


2.69 


79866 


BORA 


2.41 


55215 


FANCI 


2.68 


7072 


TIAl 


2.40 


57590 


WDFYl 


2.67 


55632 


G2E3 


2.40 


55142 


IIAUS2 


2.66 


2213 


FCGR2B 


2.40 


23047 


PDS5B 


2.66 


3987 


LIMSl 


2.39 


5373 


PMM2 


2.66 


829 


CAPZAl 


2.39 


11065 


UBE2C 


2.66 


26973 


CHORDCl 


2.38 


23085 


ERCl 


2.66 


435 


ASL 


2.38 


389197 


C4orf50 


2.65 


29979 


UBQLNl 


2.38 


11260 


XPOT 


2.65 


8548 


BLZFl 


2.37 


29980 


DONSON 


2.65 


9694 


TTC35 


2.37 


64399 


HHIP 


2.64 


55055 


ZWILCH 


2.36 


6453 


ITSNl 


2.63 


4481 


MSRl 


2.36 


29108 


PYCARD 


2.63 


10213 


PSMD14 


2.35 


9877 


ZC3H11A 


2.62 


9966 


TNFSF15 


2.35 
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Table S3. Cont. 



NCBI gene ID 


Gene Symbol 


Fold Change 


NCBI gene ID 


Gene Symbol 


Fold Change 


81624 


DIAPH3 


2.62 


51582 


AZINl 


2.35 


79723 


SUV39H2 


2.61 


54843 


SYTL2 


2.34 


55789 


DEPDCIB 


2.61 


9039 


UBA3 


2.33 


10097 


ACTR2 


2.59 


933 


CD22 


2.33 


23036 


ZNF292 


2.58 


5685 


PSMA4 


2.33 


22936 


ELL2 


2.57 


9885 


0SBPL2 


2.33 


8477 


GPR65 


2.57 


9262 


STK17B 


2.33 


23397 


NCAPH 


2.57 


56942 


C16orf61 


2.32 


3015 


H2AFZ 


2.54 


10767 


HBSIL 


2.32 


55749 


CCARl 


2.53 


87178 


PNPTl 


2.32 


25937 


WWTRl 


2.52 


6303 


SATl 


2.32 


360023 


ZBTB41 


2.51 


7316 


UBC 


2.32 


5080 


PAX6 


2.51 


4205 


MEF2A 


2.32 


4193 


MDM2 


2.51 


85465 


EPTl 


2.31 


24137 


KIF4A 


2.51 


84640 


USP38 


2.31 


9212 


AURKB 


2.51 


5810 


RADl 


2.30 


168850 


ZNF800 


2.50 


64397 


ZFP106 


2.29 


55109 


AGGFl 


2.49 


5706 


PSMC6 


2.29 


23185 


LARP4B 


2.49 


22948 


CCT5 


2.29 


51571 


FAM49B 


2.49 


10672 


GNA13 


2.29 


51077 


FCFl 


2.49 


339344 


MYPOP 


2.28 


23167 


EFR3A 


2.49 


7292 


TNFSF4 


2.28 


23468 


CBX5 


2.48 


57103 


C12orf5 


2.28 


5396 


PRRXl 


2.48 


388403 


YPEL2 


2.28 


10096 


ACTR3 


2.47 


54876 


DCAF16 


2.27 


10308 


ZNF267 


2.47 


113235 


SLC46A1 


2.27 


6782 


HSPA13 


2.47 


11177 


BAZIA 


2.27 


3832 


KIFll 


2.47 


339175 


METTL2A 


2.26 


917 


CD3G 


2.47 


26586 


CKAP2 


2.26 


80821 


DDHDl 


2.46 


55785 


FGD6 


2.26 


52 


ACPI 


2.46 


24145 


PANXl 


2.25 


4179 


CD46 


2.46 


253461 


ZBTB38 


2.25 


10499 


NC0A2 


2.44 


23232 


TBC1D12 


2.25 


60558 


GUFl 


2.44 


995 


CDC25C 


2.25 


55676 


SLC30A6 


2.43 


55974 


SLC50A1 


2.25 


6646 


SOATl 


2.43 


472 


ATM 


2.25 


5440 


P0LR2K 


2.43 


23008 


KLHDCIO 


2.24 


84955 


NUDCDl 


2.42 


10024 


TROAP 


2.24 


54739 


XAFl 


2.42 


9521 


EEFlEl 


2.24 


84295 


PHF6 


2.23 


7402 


UTRN 


2.09 


7295 


TXN 


2.23 


55589 


BMP2K 


2.08 


2710 


GK 


2.23 


158747 


M0SPD2 


2.08 


10905 


MAN1A2 


2.22 


56886 


UGGTl 


2.07 


6780 


STAUl 


2.22 


203100 


HTRA4 


2.07 
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Table S3. Cont. 



NCBIgenelD Gene Symbol Fold Change NCBI gene ID Gene Symbol Fold Change 



10282 


BETl 


2.22 


55279 


ZNF654 


2.07 


134430 


WDR36 


2.21 


54499 


TMCOl 


2.07 


4299 


AFFl 


2.21 


81930 


KIF18A 


2.07 


6747 


SSR3 


2.21 


142686 


ASB14 


2.06 


7334 


UBE2N 


2.21 


55209 


SETD5 


2.06 


5965 


RECQL 


2.21 


9736 


USP34 


2.04 


4605 


MYBL2 


2.2 


116285 


ACSMl 


2.04 


6093 


ROCKl 


2.19 


2201 


FBN2 


2.04 


161725 


0TUD7A 


2.19 


963 


CD53 


2.04 


23518 


R3HDM1 


2.18 


55159 


RFWD3 


2.03 


2239 


GPC4 


2.18 


9871 


SEC24D 


2.03 


28977 


MRPL42 


2.18 


9887 


SMG7 


2.02 


64859 


0BFC2A 


2.18 


23376 


UFLl 


2.02 


3845 


KRAS 


2.18 


79646 


PANK3 


2.01 


51388 


NIP7 


2.18 


50613 


UBQLN3 


2.00 


7586 


ZKSCANl 


2.18 


201595 


STT3B 


2.00 


10762 


NUP50 


2.17 


59345 


GNB4 


1.99 


7328 


UBE2H 


2.17 


5876 


RABGGTB 


1.99 


10730 


YMEILI 


2.17 


79820 


CATSPERB 


1.99 


23093 


TTLL5 


2.17 


6637 


SNRPG 


1.99 


6790 


AURKA 


2.17 


51330 


TNFRSF12A 


1.99 


22889 


KIAA0907 


2.17 


9928 


KIF14 


1.99 


10875 


FGL2 


2.17 


286097 


EFHA2 


1.98 


23161 


SNX13 


2.17 


9131 


AIFMl 


1.98 


9169 


SCAFll 


2.16 


488 


ATP2A2 


1.98 


1788 


DNMT3A 


2.15 


23042 


PDXDCl 


1.98 


9088 


PKMYTl 


2.15 


7114 


TMSB4X 


1.98 


23033 


DOPEY 1 


2.13 


9123 


SLC16A3 


1.98 


89882 


TPD52L3 


2.13 


54454 


ATAD2B 


1.97 


6556 


SLCllAl 


2.13 


23143 


LRCHl 


1.97 


64216 


TFB2M 


2.13 


4212 


MEIS2 


1.97 


3071 


NCKAPIL 


2.13 


1457 


CSNK2A1 


1.97 


51068 


NMD3 


2.13 


80012 


PHC3 


1.97 


509 


ATP5C1 


2.13 


128497 


SPATA25 


1.96 


953 


ENTPDl 


2.13 


186 


AGTR2 


1.96 


51105 


PHF20L1 


2.13 


53981 


CPSF2 


1.96 


5062 


PAK2 


2.13 


56996 


SLC12A9 


1.96 


9205 


ZMYM5 


2.12 


1584 


CYPllBl 


1.96 


55157 


DARS2 


2.12 


133619 


PRRCl 


1.96 


8520 


HATl 


2.11 


4288 


MKI67 


1.96 


79739 


TTLL7 


2.11 


9014 


TAFIB 


1.96 


9495 


AKAP5 


2.10 


55858 


TMEM165 


1.96 


3181 


HNRNPA2B 


2.10 


2212 


FCGR2A 


1.96 
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Table S3. Cont. 



NCBIgenelD Gene Symbol Fold Change NCBI gene ID Gene Symbol Fold Change 



389898 


UBE2NL 


2.10 


10075 




HUWEl 


1.96 


29850 


TRPM5 


2.10 


220988 




HNRNPA3 


1.96 


3070 


HELLS 


2.10 


80146 




UXSl 


1.95 


331 


XIAP 


2.09 


122011 




CSNKIAIL 


1.95 


55751 


TMEM184C 


2.09 


150468 




CKAP2L 


1.95 


2146 


EZH2 


2.09 


84624 




FNDCl 


1.95 


26057 


ANKRD17 


1.95 


7332 




UBE2L3 


1.92 


128061 


Clorfl31 


1.95 


3336 




HSPEl 


1.92 


64090 


GAL3ST2 


1.94 


54800 




KLHL24 


1.92 


130507 


UBR3 


1.93 


2290 




FOXGl 


1.91 


2298 


F0XD4 


1.93 


50848 




FUR 


1.91 


123169 


LEOl 


1.93 


10627 




MYL12A 


1.91 


57187 


TH0C2 


1.93 


5074 




PAWR 


1.91 


148789 


B3GALNT2 


L93 


6476 




SI 


1.91 


58508 


MLL3 


1.92 


1009 




CDHll 


1.90 


5701 


PSMC2 


1.92 


29066 




ZC3H7A 


1.90 


148066 


ZNRF4 


1.92 


51319 




RSRCl 


1.90 


6670 


SP3 


1.92 










Table S4. Result of ROC curve analysis on our miRNA array data. ROC analysis was done 


to validate the diagnostic value of the miRNA in the miRNA-regulated PINs. AUC: area 


under (ROC) curve; */7-value < 0.05; 


*H=*^.value< 0.001. 








miRBase Accession 


miRNA name 


AUC 


/)-value 






MIMAT0002856 


hsa-miR-520d-3p 


0.49 


0.549112 






MIMAT0000265 


hsa-miR-204-5p 


0.98 


6.47 X 10"'° *** 






MIMAT0000272 


hsa-niiR-215 


0.21 


0.999782 






MIMAT0000271 


hsa-miR-214-3p 


0.68 


0.010387 * 






MIMAT0002820 


hsa-miR-497-5p 


0.99 


2.75 X 10"'° *** 






MIMAT0000076 


hsa-miR-2 1 -5p 


0.78 


0.000184 






MIMAT0000738 


hsa-miR-383 


0.60 


0.106284 






MIMAT0000423 


hsa-miR-125b-5p 


0.99 


2.48 X 10"'° *** 






MIMAT0000064 




hsa-let-7c 


0.93 


3.79 X 10"^ *** 






MIMAT0000089 


hsa-miR-31-5p 


0.80 


8.63 X 10"^ 






MIMAT0000077 


hsa-miR-22-3p 


0.27 


0.99749 






MIMAT0000098 


hsa-miR-100-5p 


0.98 


5.55 X 10"'° *** 






MIMAT0000097 


hsa-niiR-99a-5p 


0.99 


2.55 X 10"'° *** 






MIMAT0000443 


hsa-mlR- 1 25a-5p 


0.31 


0.990694 






MIMAT0002819 


hsa-miR-193b-3p 


0.41 


0.86128 






MIMAT0000250 


hsa-miR-139-5p 


0.99 


2.42 X 10"'° *** 






MIMAT0000437 


hsa-miR-145-5p 


0.96 


3.14 X 10"^ *** 






MIMAT0000421 


hsa-miR-122-5p 


0.48 


0.597483 
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Table S5. Summary of constructed miRNA-regulated networks. LO gene: genes connected 
directly to the iniRNA {i.e., direct target of miRNA); LI gene: genes not connected directly 
to the miRNA. 



miRBase Accession 


miR name 


Total gene count 


LO count 


LI count 


MIMAT0002819 


hsa-miR-193b-3p 


16 


1 


15 


MIMAT0000250 


hsa-miR-139-5p 


28 


10 


18 


MIMAT0000437 


hsa-miR-145-5p 


86 


22 


64 


MIMAT0000423 


hsa-miR-125b-5p 


211 


16 


195 


MIMAT0000443 


hsa-miR-125a-5p 


206 


14 


192 


MIMAT0000097 


hsa-miR-99a-5p 


14 


1 


13 


MIMAT0000265 


hsa-miR-204-5p 


64 


18 


46 


MIMAT0000076 


hsa-niiR-21-5p 


91 


16 


75 


MIMAT0000064 


hsa-let-7c 


96 


20 


76 


MIMAT0000421 


hsa-miR-122-5p 


5 


3 


2 


MIMAT0000098 


hsa-miR-100-5p 


14 


1 


13 


MIMAT0000272 


hsa-miR-215 


3 


3 


0 


MIMAT0000271 


hsa-miR-214-3p 


14 


8 


6 


MIMAT0000738 


hsa-miR-383 


34 


3 


31 


MIMAT0002856 


hsa-miR-520d-3p 


146 


23 


123 


MIMAT0000077 


hsa-miR-22-3p 


46 


11 


35 


MIMAT0002820 


hsa-miR-497-5p 


267 


32 


235 


MIMAT0000089 


hsa-miR-31-5p 


34 


3 


31 
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Table S6. Specific enriched GO terms of each miRNA-reguIated PINs. Genes annotated with the specific GO term in the PIN were also listed 
in this table. Adj.jc-value: multiple-test adjusted p-value calculated by the method described in the work of Benjamini and Yekutieli [48]. 

MIMAT0002856(hsa-miR-520d-3p) 



GO term 



Genes 



Adj.^-value 



^169, Transmembrane SH3KBP1, HDAC2, RET, ABIl, LYN, GRB2, SORBSl, CLTC, CLTA, CDC42, CASP9, RAFl, SRC, AP2A1, 

atein tyrosine kinase AP2B1, MAPK3, ARHGEF7, PRKCA, RPS6, PRKAR2B, MAPKl, ARHGEF6, CDKl, SH3GL2, EIF4G1, HDACl, 2.77 x If 



GO:004801 1, Nerve growth factor 
receptor signaling patliway 
GO:0007173, Epidermal growth 
factor receptor signaling path wa.y^ 



HDAC2, GRB2, CLTC, CLTA, CASP9, RAFl, SRC, AP2A1, AP2B1, MAPK3, ARHGEF7, PRKCA, PRKAR2B, 
MAPKl, ARHGEF6, CDKl, SH3GL2, HDACl, ECT2, CASP3, PRKACA, PRKAR2A, SHCl, RACl 
SH3JCBP1, GRB2, CLTC, CLTA, CDC42, CASP9, RAFl, SRC, AP2A1, AP2B1, MAPK3, ARHGEF7, PRKCA, 
PmAR2R.MA£KL CD.Kl„SmGO...BRKACA..ERKAR? A SHCl 



5.09 X 10" 



2.96 X 10 



GO:0043067, Regulation of 
programmed cell death 

GO:0042058, Regulation of 
epidermal growth factor 
llfig^^tg ^ i ^nalin g pathway 
GO:0008543, Fibroblast growth 
factor receptor signaling pathway 
GO:0043068, Positive regulation 
of programmed cell deat^ 
GO:0010942, Positive 
regulation of cell death 



HDAC2, STK17B, ESRl, ABLl, LYN, TP53, GABRB3, PAK2, LCK, CASP9, RAFl, PLKl, ARHGEF7, PRKCA, 
RPS6, SH3RF1, MAPKl, IFT57, ARHGAPIO, ARHGEF6, CDKl, APAFl, HDACl, ECT2, CASP3, SOXIO, EP300, 
ARAF, TFAP2A, ADRB2, HCK, KLHL20, CASP8, HIPl, RACl 



SH3KBP1, ESRl, GRB2, CLTC, CLTA, CDC42, AP2A1, AP2B1, ARHGEF7, SH3GL2, SHCl 




4.76 X 10"'' 



2.36 X 10"'^ 



GRB2, CASP9, RAFl, SRC, MAPK3, PRKCA, PRKAR2B, MAPKl, CDKl, MKNKl, PRKACA, PRKAR2A, SHCl 2.39 x 10" 



STK17B, ABLl, LYN, TP53, LCK, CASP9, ARHGEF7, PRKCA, RPS6, SH3RF1, MAPKl, ARHGEF6, APAFl, 
ECT2, CASP3, EP300, TFAP2A, ADRB2, CASP8, HIPl, RACl 

STK17B, ABLl, LYN, TP53, LCK, CASP9, ARHGEF7, PRKCA, RPS6, SH3RF1, MAPKl, ARHGEF6, APAFl, 
ECT2, CASP3, EP300, TFAP2A, ADRB2, CASP8, HIPl, RACl 



4.49 X 10" 



STK17B, ABLl, TP53. LCK, CASP9, ARHGEF7, PRKCA, SH3RF1, MAPKl, ARHGEF6, APAFl, ECT2, CASP3, 




GO :00 12502, Induction of 
programmed cell death 



STK17B, ABLl, TP53, LCK, CASP9, ARHGEF7, PRKCA, SH3RF1, MAPKl, ARHGEF6, APAFl, ECT2, CASP3, 
EP300, CASP8, HIPl, RACl 



7.42 X 10"" 
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Table S6. Cont. 



MIMAT0002856(hsa-miR-520d-3p) 


GO term 


Genes 


Adj.p-value 


KT ^ 

00:0042059, Negative regulation 






of epidermal growth factor 


SH3KBP1, ORB2, CLTC, CLTA, CDC42, AP2A1, AP2B1, ARHOEF7, SH30L2 


8.63 X 10"" 








00:0015630, Microtubule 

cytoskeleton 


STMNl, SORBSl, SMAD4, CLTC, CDC42, LCK, RACOAPl, PLKl, PRKAR2B, YESl, MAPKl, IFT57, CDKl, 

ECT2, PRKACA, RBI, EP300, CCNBl, CHAFIB, TFAP2A, CASP8, PRKAR2A 


4.75 X 10"'" 


00:0060548, Negative regulation 


HDAC2, ESRl, TP53, SMAD4, RAFl, PLKl, PRKCA, RPS6, SH3RF1, CDKl, HDACl, CASP3, SOXIO, ARAF, 


6.27 X 10"* 


00:0008286, Insulin receptor 

signaling pathway 


GRB2, SORBSl, RAFl, MAPK3, RPS6, MAPKl, CDKl, EIF4G1, EIF4B, SHCl 


2.15 X 10"' 


00:0043069, Negative regulation 


HDAC2, ESRl, TP53, RAFl, PLKl, PRKCA, RPS6, SH3RF1, CDKl, HDACl, CASP3, SOXIO, ARAF, TFAP2A, 


3.13 X 10"' 


iof nroerammed cell death 


HCK.KLHL20_ 


00:0008284, Positive regulation 
of cell proliferation 


HDAC2, ESRl, LYN, CDC42, E2F1, PRKCA, MAPKl, CDKl, RHOO, HDACl, NCKl, SOXIO, CCNBl, ADRB2, 

HCK, SHCl 


8.34 X 10"' 


00:0051988, Regulation of 






littachment of spindle 


CDC42, RACOAPl, ECT2, ccNBl.^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^B 




microtubules to MtttflStftiBii^^H 






00:0008629, Induction of 
apoptosis by intracellular signals 


ABLl, TP53, CASP9, APAFl, CASP3, EP300, CASP8 


5.83 X 10"^ 


MIMAT0002820(hsa-miR-497-5p) 


GO term 


Genes 


Adj./>-value 




243067, Regulation of 
ed cell death 



00:0043068, Positive regulation 
of programmed cell death 

00:0010942, Positive regulation 



2.67 X 10 




ESRl, MENl, ABLl, HIPK3, PPAROCIA, SIAin^H3RFl, PAK2, LCK, MEDl, PPARO, CBX4, ARHOEF7, 
YWHAB, RXRA, ACVRl, MAPKl, CASP3, CASP6, AR, PTPRF, MDM2, BRCAl, MLHl, RAB27A, PIAS4, FAFl, 
RACl, VHL, SKI, NR4A1, LYN, TP53, PSMC2, GATAl, GATA6, 0ATA3, RAFl, CDKNIB, PLKl, PSMDll, 
HOXA13, RPS6, ESR2, ARHOAPIO, ARHOEF6, SMAD3, SKIL, RYR2, PSENl, HCK, TRAF2 
MENl, ABLl, SIAHl, SH3RF1, LCK, PPARO, ARHOEF7, YWHAB, RXRA, MAPKl, CASP3, CASP6, PTPRF, 
BRCAl, MLHl, RAB27A, PIAS4, FAFl, RACl, NR4A1, LYN, TP53, OATA6, CDKNIB, H0XA13, RPS6, ESR2, 1.74 x 10"" 
ARHGEF6, SMAD3, RYR2, PSENl, TRAF2 

MENl, ABLl, SIAHl, SH3RF1, LCK, PPARO, ARHOEF7, YWHAB, RXRA, MAPKl, CASP3, CASP6, PTPRF, 
MLH1,£A^2A. PIAS4, FAFl, RACl, NR4A1, LYN, TP53, OATA6. CaKJiJ:m--HQX.All-£LES6. ESR2. 
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Table S6. Cont. 



MIMAT0002820(hsa-miR-497-5p) 



GO term 



Genes 



Adj.p-value 

jO:0008285, Negative regulation MENl, MEDl, PPARG, RXRA, CASP3, AR, PTPRF, VDR, VHL, SKI, LYN, TP53, tOBl, GATAl, GATA3, RAFl, 3 ^ , „-i4 

F4A, CDKNIB, BRD7, MED25, ESR2, ABIl, SMADl, SMAD2, SMAD3, SMAD4, S0X7 




GO:00 15629, Actin cytoskeleton 
}O:00069T7r 




ABLl, SORBSl, FLNA, SEPT7, ANLN, MACFl, HAPl, SH3PXD2A, 1QGAP2, BRCAl, ACTCl, ACTAl, MYL2, .3 
MYLK, SORBS2,ARPC4,ARPC5,ACTR2,ACTR3,ARPClB, WASFl, WASF2, HCK .5x0 
ABLl, SH3RF1, LCK, PPARG, ARHGEF7, YWHAB, MAPKl, CASP3, CASP6, BRCAl, MLHl, RAB27A, RACl, 



GO:00 12502, Induction of 
programmed cell death 
GO:0007178, Transmembrane 
receptor protein serine/threonine 
Mnas e signaling pathway 
GO:0007169, Transmembrane 
receptor protein tyrosine kinase 
signaling pathway 

390092, Regulation of 
transmembrane receptor protein 
serine/threonine kinase signaling 
■Dathmav, 

GO:0030509, BMP signaling 

pathway 



TP53- CDKN 1B ARHGEF6, SMAD3, RYR2, PSENl, TRAF2 
ABLl, SH3RF1, LCK, PPARG, ARHGEF7, YWHAB, MAPKl, CASP3, CASP6, BRCAl, MLHl, RAB27A, RACl, 
NR4A1, TP53, CDKNIB, ARHGEF6, SMAD3, RYR2, PSENl, TRAF2 




9.6' 



1.06 X 10" 



ACVRl, SMURF2, SKI, GDF6, BMP6, ZNF8, GATA4, HNF4A, SMADl, SMAD2, SMAD3, SMAD4, SMAD5, ^^_io 
RYR2 ^-^^"^^ 



SORBSl, CDC42, SRC, MAPK3, ARHGEF7, YWHAB, MAPKl, SH3GL2, CASP3, MDM2, EIF4G1, RACl, 10 
SH3KBP1, NR4A1, LYN, GRB2, RAFl, CDKNIB, RPS6, ABIl, ARHGEF6, MKNKl, PSENl, EIF4B 




ACVRl, SMURF2, SKI, GDF6, BMP6, ZNF8, SMADl, SMAD4, SMAD5, RYR2 



3.25 X 10" 




GO:0007173, Epidermal growth CDC42, SRC, MAPK3, ARHGEF7, YWHAB, MAPKl, SH3GL2, MDM2, SH3KBP1, NR4A1, GRB2, RAFl, 9 22 x 10"' 
factor receptor signaling pathway CDKN 1 B 
Lis 



GO:0043069, Negative regulation ESRl, HIPK3, PPARGCIA, SH3RF1, MEDl, CBX4, ACVRl, CASP3, AR, MDM2, VHL, TP53, GATAl, GATA6, g 
of program med cell death GATA3 , RAF 1 , CDKN 1 B, PLKl , RPS6, SMAD3 , PSEN 1 , HCK 

PmUBl 1, Nerve growth factor^ 
receptor signaling pathway 



GO:0032956, Regulation of actin 
cytoskeleton organization 



ABLl, LRPl, ARPC4, ARPC5, ACTR3, ARPCIB, SMAD3, NCKl, S0RBS3, HCK, LIMKl 



6.42 X 10" 
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MIMAT0002820(hsa-iniR-497-5p) 



GO term 



Genes 



Adj./j-value 



GO:0008543, Fibroblast growth 




SRC, MAPK3. YWHAB. MAPKl. MDM2. NR4A1. GRB2. RAFl. CDKNIB. MKNK 



factor receptor signaling pathway | 
GO:0042059, Negative regulation 

of epidermal growth factor CDC42, ARHGEF7, SH3GL2, PTPRF, SH3KBP1, GRB2, PSENl 

receptor signaling pathway 
GO:0042058, Regulation of 
epidermal growth factor receptor ESRl , CDC42, ARHGEF7, SH3GL2, PTPRF, SH3KBP1, GRB2, PSENl. 



1.11 X 10"^ 



1.17 X 10" 



ACVRl, SMURF2, SKI, SMADl, SMAD2, SMAD3, SMAD4, SMAD5 



STMNl, RIFl, SORBSl, CDC42, LCK, RACGAPl, YESl, YWHAB, MAPKl, SEPT7, KIF23, CDC16, MACFl, 



GO:0007179, Transforming 
growth factor beta receptor 
signaling pathway 
trO:00 15630, Microtubule 
cytoskeleton 
GO:0017015, Regulation of 

transforming growth factor beta MENl , SMURF2, SKI, TP53, SMAD2, SMAD3, SMAD4, SKIL 

receptor signaling pathway 
GO;0070302, Regulation of 



1.67 X 10" 



2.30 X 10" 




HSP90AB1, MEDl, PPARG, PTPRF, NR1H2, PIASl, IL36RN, HIPKl 



6.35 X 10" 



GO:0001959, Regulation of 
cytokine-mediated 

signaling pathway 

[gO:0008284, Positive regulation ESRl, CDC42, MEDl, RARA, MAPKl, AR, MDM2, NR4A1, LYN, FZRl, BMP6, GATAl, GATA4, GATA6, „ ,._5 
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MIMAT0000423(hsa-miR-125b-5p) 



GO term 



Genes 



Adj.p-value 



GO:0043067, Regulation of 
programmed cell death 



GO:0060548, Negative regulation 
of cell death 

GO:0008284, Positive regulation 
'cellprohferation 



4.56 X 10 



■24 




HMGA2, PML, PRNP, FGF2, XRCC4, BRCAl, IGFBP3, HDAC3, CTNNBl, CDS, CDKl, NKX2-5, MEF2C, PRKCI, 

CASP2, PSMA4, PSMA3, CFDPl, CAVl, FAFl, YWHAB, HIFIA, RELA, TCF7L2, TNFSF12, PSEN2, TP53, 
TOP2A, TNFRSF4, BID, MYC, JUN, OGT, CDKNIA, HOXA13, RNF7, PPP2R4, HDAC2, HDACl, SNCA, PTEN, 
NFKBIA, IFI16, N0L3, TRAF2, HSP90B1 

HMGA2, PRNP, FGF2, XRCC4, HDAC3, CTNNBl, CDKl, NKX2-5, MEF2C, PRKCI, CFDPl, HIFIA, RELA, 
TCF7L2, PSEN2, TP53, TNFRSF4, MYC, JUN, CDKNIA, RNF7, HDAC2, HDACl, SNCA, PTEN, MGMT, 8.04 x lo '^ 
NFKBIA, NOL3 , HSP90B 1 _ 
HMGA2, FGF2, XRCC4, CDC25B, CTNNBl, EGRl, AGGFl, CDKl, NKX2-5, MEF2C, PRKCI, IRSl, HIFIA, 
RELA, HCLSl, TNFSF12, ARNT, PTPRC, TNFSF4, TNFRSF4, MYC, JUN, FGFl, CDKNIA, HDAC2, HDACl, 2.20 x lO ^^ 
NOLCl,PTEN 

143069, Negative regulatio n HMGA2. PRNP. XRCC4, HDAC 3. CT^^s^RL. CDK.L. MKX2^S .MEF2C- PRt^^ CFDPl. HIFIA, REL A, TCF7L2, 
ledcell d^th, 

GO:0043068, Positive regulation HMGA2, PML, BRCAl, IGFBP3, CTNNBl, CDS, MEF2C, PRKCI, CASP2, CAVl, FAFl, YWHAB, TNFSF12, 
of program med cell death PSEN2, TP53, TOP2A, BID, JUN, OGT, CDKNIA, HOXA13, RNF7, PPP2R4, PTEN, IFI16, TRAF2 

Positiveregulation HMGA2, PML, BRCAl, IGFBP3, CTNNBl, CDS, MEF2C, PRKCI, CASP2, CAVl, FAFl, YWHAB, TNFSF12, 
PSEN2, TP53, T0P2A, BID, JUN, OGT, CDKNIA, H0XA13, RNF7, PPP2R4, PTEN, IFI16, TRAF2 
PRNP, HDAC3, CDKl, NKX2-5, MEF2C, PRKCI, CFDPl, RELA, TCF7L2, PSEN2, RNF7, HDACl, SNCA, 
NFKBLV, NOL3, HSP90B1 

85, Negative regulation SERPINFl, SRF, PML, PRNP, FGF2, CSNK2B^ IGFBP3, CTNNBl, CAVIVhMGAI, VDR, CDH5, HSFl, C0L18A1^ . 
of cell prohferation TP53, MYC, JUN, CDKNIA, PAKl , PTEN ^g||||^j^^ 



4.51 X 10-1" 




GO:0006916, Anti-apoptosis 



1.25 X 10 ' 



GO:0015630, Microtubule 
cytoskeleton 



STMNl, KIFIC, RANGAPl, CDC25B, BRCAl, HDAC3, CTNNBl, PIN4, HSPHl, RANBP9, CDKl, SPTANl, 
YWHAQ, DVLl, FKBP4, YWHAB, CCDC85B, MAPT, PSEN2, T0P2A, SPIB, MYC, OGT, APEXl, PAFAHIBI 



2.24 X 10"' 



pO:004801 1, Nerve growth factor 
receptor signaling pathway 
GO:0050678, Regulation of 
epithelial cell proliferation 



CDKl. MEF2C. PRKCI. CASP2. IRSl. YWHAB. RELA. PSEN2. HDAC2. HDACl. PTEN. ATFl. NF 



SERPINFl, PGR, FGF2, CTNNBl, AGGFl, CAVl, HIFIA, TNFSF12, ARNT, MYC, JUN, FGFl, PTEN 



3.41 X 10" 



GO:0007 169, Transmembrane ^^^^ 

FGF2, 

receptor protein tyrosme kinase , ^, 

. , HDACl 

l^naling pathway 
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MIMAT0000423(hsa-miR-125b-5p) 


GO term 


Genes 




Adj. /J- value 




00 0012502 Tndiictinn of 
programmed cell death 


PML, BRCAl, CDS, CASP2, CAVl, YWHAB, TNFSF12, PSEN2, TP53, BID, OOT, CDKNIA, RNF7, PTEN, IFI16, TRAF2 


1.81 X lO"' 


^^^035666 TRTF-denendent 








PHHP^ceptor signaling 

ipaUlWay 


ATF2, CDKl, MEF2C, FOS, RELA, JUN, ATFl, NFKBL\ 






00:0034138, Toll-like receptor 3 


ATF2, CDKl, MEF2C, FOS, RELA, JUN, ATFl, NFKBIA 




2.71 X 10"* 


^fttKf/k93, Actin filament 
capping 

00:0002756, MyD88- 


SPTB, SPTBNl, SPTANl, SPTAl, ADDl, EPB49 ^^^H 




3.43 X 


independent toll-like recegtor 
igignaling pathway 


ATF2, CDKl, MEF2C, FOS, RELA, JUN, ATFl, NFKBLV 




3.82 X 10"* ^ 


00:0034134, Toll-like receptor 2 

signaling pathway 


ATF2, CDKl, MEF2C, FOS, RELA, JUN, ATFl, NFKBIA 




5.33 X 10"* 


00:0034130, Toll-like receptor 1 


ATF2, CDKl, MEF2C, FOS, RELA, JUN, ATFl, NFKBIA 






00:0030835, Negative regulation 
of actin filament depolymerization 


SPTB, SPTBNl, SPTANl, SPTAl, ADDl, EPB49 




5.58 X 10"* 


trO:0015629, cytoskeleton^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^BHIHHIi| 


00:0002755, MyD88-dependent 








toll-like receptor signaling 


ATF2, CDKl, MEF2C, FOS, RELA, JUN, ATFl, NFKBIA 




7.66 X 10"* 


pathway 








EO:0034142, Toll-i'ike receptor 4 
[signaling pathway | 


,4ffiaMQ^LMEE2g. FQS. Rela, jun, atfi,^^^^^^ 






00:0050679, Positive regulation 
of epithelial cell proliferation 


FOF2, CTNNBl, AOOFl, HIFIA, TNFSF12, ARNT, MYC, JUN, FOFl 


1.14 X 10"' 
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MIMAT0000076(hsa-miR-21-5p) 


GO term 


Genes 


Adj./>-value 


GO:0030834, Regulation of actin 
filament depolymerization 


SPTB. SPTBNl, SPTANl, SPTAl, APr^' rpp/iq^^^^^^^^^^^^^^^^^^^^^^^^^^^B 


Ll.27 X 10"^ ^ 


GO:0030837, Negative regulation 
of actin filament polymerization 


SPTB, SPTBNl, SPTANl, SPTAl, ADDl, EPB49 


2.71 X 10"' 


pO:0002224, Toll-like receptor 


ATF7 rr)K1 MFF9r FDS RFT A HTN ATF1 NFKRTA 


2.96 X 10-5 


GO:0008629, Induction of 
apoptosis by intracellular signals 


PML, BRCAl, YWHAB, TP53, BID, CDKNIA, RNF7, IFI16 


3.52 X 10"' 


00:0043067, Regulation of 
programmed cell death 


SPRY2, TP53, ADAMTSL4, ETSl, TDGFl, RAFl, H0XA5, H0XA13, MSXl, MSX2, NKX2-5, CBL, INHBB, 

COL4A3, ACVRIC, TRAF2 

^^^^^^^^^^Ull, TDGFl, RAFl, CBL ^^^^^^^^^^^^^^^^^^^^ 


1.92 X 10"' 1 


MIMAT0000250(hsa-miR-139-5p) 


GO term 


Genes 


Adj.p-value 



HPPBIi83, Negative regulation 
pf transforming growth factor beta 
receptor signaling pathway by LTBP 1 , FBN 1 , FBN2 
extracellular sequestering of 




MIMAT0000089(hsa-miR-31-5p) 



GO term 



Genes 




0:0007187, G-protein signaling, 

oupled to cychc nucleotide GNA12, GNA13, DRDS, MTNR1A, S1PR3, TSHR, S1PR4 

second messenger 
GO:0019935, Cyclic-nucleotide- 
mediated signaling 
^0:0007188, G-protein signaling, 

coupled to cAMP nucleotide GNA12, GNA13, DRD5, S1PR3, TSHR, S1PR4 
second messenger — 



Adj. j?-value 



GNA12, GNA13, DRD5, MTNRIA, S1PR3, TSHR, S1PR4 



1.60 X 10"^ 
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MIMAT0000089(hsa-miR-31-5p) 



GO term 



Genes 



Adj. j?-value 



GO:004801 1, Nerve growth factor 

receptor signaling pathway 
GO:00 19933, CAMP-mediated 
Igjgnaling _ 
GO:0043067, Regulation of 
programmed cell death 
jO;0003376, Sphingosine-l 



ARHGEFl, PRKCD, ARHGEF12, PRKACA, PRKCE, MCF2, ARHGEFl 1 
GNA12, GNA13, DRD5, S1PR3, TSHR, S1PR4 
PTK2B, PRKCD, TGFBRl, ARHGEF12, CTNNBl, PRKCE, TIAl, MCF2, FASTK, ARHGEFll, F2R 



GO:0007169, Transmembrane 

receptor protein tyrosine kinase PTK2B, ARHGEFl, PRKCD, ARHGEF12, PRKACA, PRKCE, MCF2, ARHGEFl 1 
signaling pathway 
||rO:0043068, Positive regulation 
of programmed cell death 
GO:00 10942, Positive regulation 
of cell death 

K):0006917, Induction of 
apoptosis 



TGFBRl, ARHGEF12, CTNNBl, PRKCE, TIAl, MCF2, FASTK, ARHGEFl 1 
TGFBRl, ARHGEF12, CTNNBl, PRKCE, TIAl, MCF2, FASTK, ARHGEFll 
TGFBRl, ARHGEF12, PRKCE, TIAl, MCF2, FASTK, ARHGEFl 1 



MIMAT0000437(hsa-miR-145-5p) 



6.93 X 10 




1.25 X 10" 




5.20 X 10" 




8.78 X 10" 



9.32 X 10" 




GO term 



Genes 



GO:0007178, Transmembrane 

receptor protein serine/threonine BMP6, ZNF8, ACVRl, SMADl, SMAD4, RYR2, SMAD5, SMURF2, GDF6 

kinase signaling pathway 

GO:0090092, Regulation of -----------^ 



transmembrane receptor protein 
serine/threonine kinase signaline 



MENl, TP53, BMP6, H0XA13, ACVRl, SMAD4, SULFl, SMURF2, GDF6 



Adj./?-value 




2.56 X 10" 
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MIMAT0000064(hsa-let-7c) 



GO term 



Genes 



Adj./>-value 



D43067, Regulation of 
programmed cell death 
GO:0043069, Negative regulation 
of programmed cell death 
IjO:0060548, Negative regulation 

GO:0015630, Microtubule 
cytoskeleton 

GO:0043067, Regulation of 
progra 



00:0060548, Negative regulation 
of cell death 




RRM2B, ACV Rl, TP53, RASAl , TOFBRl, PSM A3, BIRC5, ACTN2, H0XA13, IRS2, FASTK, VAVl, PSMB6, 

RRM2B, ACVRl, TP53, RASAl, TGFBRl, BIRC5, IRS2, BCL2, CDKl, HDACl, SOXIO, AKTl, AURKB 

MJyOB, ACVRl, TP53, R ASAl, TGFBRl, BIRC5, IRS2, BCL2, CDKl, HDACl, SOXIO, AKTl, AURKB 

INCENP, SNTB2, SEPTl, TACCl, BIRC5, RACGAPl, PIN4, CDCA8, CDKl, PHFl, AKTl, AURKB, NINL, 
CCDC85B 

HMGA2, PML, PRNP, FGF2, XRCC4, BRCAl, IGFBP3, HDAC3, CTNNBl, CD5, CDKl, NKX2-5, MEF2C, PRKCI, 
CASP2, PSMA4, PSMA3, CFDPl, CAVl, FAFl, YWHAB, HIFIA, RELA, TCF7L2, TNFSF12, PSEN2, TP53, 
T0P2A, TNFRSF4, BID, MYC, JUN, OGT, CDKNIA, RNF7, PPP2R4, HDAC2, HDACl, SNCA, PTEN, NFKBIA, 
IFI16, NOL3, TRAF2, HSP90B1 

HMGA2, PRNP, FGF2, XRCC4, HDAC3, CTNNBl, CDKl, NKX2-5, MEF2C, PRKCI, CFDPl, HIFIA, RELA, 
TCF7L2, PSEN2, TP53, TNFRSF4, MYC, JUN, CDKNIA, RNF7, HDAC2, HDACl, SNCA, PTEN, MGMT, 
NFKBIA, NOL3, HSP90B1 

HMGA2, FGF2, XRCC4, CDC25B, CTNNBl, EGRl, AGGFl, CDKl, NKX2-5, MEF2C, PRKCI, IRSl, HIFIA, 
RELA, HCLSl, TNFSF12, ARNT, PTPRC, TNFSF4, TNFRSF4, MYC, JUN, FGFl, CDKNIA, HDAC2, HDACl, 



6.58 X 10" 



8.92 X 10" 



5.15 X 10" 



1.62 X 10 



-23 



4.53 X 10" 



1.23 X 10" 



MIMAT0000443(hsa-miR-125a-5p) 



GO term 



Genes 



Adj./>-value 



3^ 



00:0043069, Negative regulation 
Ijif _programmed cell death 
00:0043068, Positive regulation 
of programmed cell death 

^2, Positive regulation 



HMGA2, PRNP, XRCC4, HDAC3, CTNNBl, CDKl, NKX2-5, MEF2C, PRKCI, CFDPl, HIFIA, RELA, TCF7L2, 
PSEN2, TP53, TNFRSF4, MYC, JUN, CDKNIA, RNF7, HDAC2, HDACl, SNCA, PTEN, N FKBIA, NOL3. HSP90B1- 
HMGA2, PML, BRCAl, IGFBP3, CTNNBl, CD5, MEF2C, PRKCI, CASP2, CAVl, FAFl, YWHAB, TNFSF12, 
PSEN2, TP53, T0P2A, BID, JUN, OGT, CDKNIA, RNF7, PPP2R4, PTEN, IFI16, TRAF2 

iGA2. PML. BRCAl, IGFBP3. CTNNBl, CD5. MEF2C. PRKCI. CASP2. CAVl. FAFl. YWHAB. TNFSFl 



2.02 X 10 



2.80 X 10" 



. ■ PRNP, HDAC3, CDKl, NKX2-5, MEF2C, PRKCI, CFDPl, RELA, TCF7L2, PSEN2, RNF7, HDACl, SNCA, „ 

00:0006916, Anti-apoptosis ^FKBLV, NOL3, HSP90B1 

285. Negative regulalicm-^SLEINFL SRF. PML. PRNP. FGF2. CSNK2B. IGFBP3. CTNNBL CAVL HMGAl. VDR. CDH5. HSFl. COL18AL 
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MIMAT0000443(hsa-miR-125a-5p) 



GO term 



Genes 



Adj./>-value 



GO:00 15630, Microtubule 

cytoskeleton 



STMNl, KIFIC, RANGAPl, CDC25B, BRCAl, HDAC3, CTNNBl, PIN4, HSPHl, RANBP9, CDKl, SPTANl, ^ , 

YWHAQ, DVLl, FKBP4, YWHAB, CCDC85B, MAPT, PSEN2, T0P2A, SPIB, MYC, OGT, APEXl, PAFAHIBI 




00:0050678, Regulation of 
epithelial cell proliferation 
fciO:0007169, Transmembrane 
peceptor protein tyrosine kinase 
signaling pathway 
00:0006917, Induction of 

apoptosis 

^^^12502, Induction of 
pBgrammed cell death ^ 



SERPINFl, POR, FOF2, CTNNBl, AOOFl, CAVl, HIFIA, TNFSF12, ARNT, MYC, JUN, FOFl, PTEN 

FOF2, HDAC3, CDKl, MEF2C, PRKCl, CASP2, IRSl, FIBP, PTPNl, YWHAB, RELA, PSEN2, FOFl, HDAC2, 
HDACl, PTEN, ATFl, NFKBIA, E1F4EBP1 

PML, BRCAl, CD5, CASP2, CAVl, YWHAB, TNFSF12, PSEN2, TP53, BID, OOT, CDKNIA, RNF7, PTEN, 1F116, 

TRAF2 

PML, BRCAl, CD5, CASP2, CAVl, YWHAB, TNFSF12, PSEN2, TP53, BID, OOT, CDKNIA, RNF7, PTEN, IFI16, 




00:0035666, TRIF-dependent 
toll-like receptor signaling 

pathway 

00:0034138, Toll-like receptor 3 
^naling pathway ^^^^^^^^^ 
00:0051693, Actin filament 

capping 

00:0002756, MyD88- 
independent toll-like receptor 
signaling pathway 
00:0015629, Actin cytoskeleton 



ATF2, CDKl, MEF2C, FOS, RELA, JUN, ATFl, NFKBIA 



2.01 X 10"^ 



ATF2, CDKl, MEF2C, FOS, RELA, JUN, ATFl, NFKBIA 



SPTB, SPTBNl, SPTANl, SPTAl, ADDl, EPB49 




2.97 X 10" 



ATF2, CDKl, MEF2C, FOS, RELA, JUN, ATFl, NFKBIA 

WAS, CDHl, BRCAl, SPTB, SPTBNl, SPTANl, CTDPl, STXIA, SPTAl, PAKl, SNCA, ADDl, EPB41, EPB49 




4.25 X 10" 




00:0034130, Toll-like receptor 1 
signaling pathway 



ATF2, CDKl, MEF2C, FOS, RELA, JUN, ATFl, NFKBIA 



4.43 X 10"^ 
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MIMAT0000443(hsa-miR-125a-5p) 


GO term 


Genes 


Adj.p-value 


GO:0002755, MyD88-dependent toll-like receptor signaling pathway 


ATF2, CDKl, MEF2C, FOS, RELA, JUN, ATFl, NFKBIA 


6.41 X 10"'' 


GO:0050679, Positive regulation of epithelial cell pr^^^^^^^^H 


1 FGF2, CTNNBl, AGGFl, HIFIA, TNFSF12, ARNT, MYC, i^^^l^^ 


^^^^^9.31 X loH 

9.31 X 10"* 


GO:0034142, Toll-like receptor 4 signaling pathway 


ATF2, CDKl, MEF2C, FOS, RELA, JUN, ATFl, NFKBIA 


|gO:0030834, Regulation of actin filament depolymerization 


SPTB, SPTBNl, SPTANl, SPTAl, ADDl, EPB49 


i.09^vm 


GO:0030837, Negative regulation of actin filament polymerization 


SPTB, SPTBNl, SPTANl, SPTAl, ADDl, EPB49 


2.37 X 10"^ 


GO:0002224, Toll-like receptor signaling pathway tfl^H 


^ ATF2, CDKl, MEF2C, FOS, RELA, JUN, ATFl, NFKBIA JHHHi 


■■■■111^2.48 X lOjig 

2.92 X 10"^ 


GO:0008629, Induction of apoptosis by intracellular signals 


PML, BRCAl, YWHAB, TP53, BID, CDKNIA, RNF7, IFI16 




GO:0002221, Pattern recognition receptor signaling pathway 


ATF2, CDKl, MEF2C, FOS, RELA, JUN, ATFl, NFKBIA 


9.15 X 10"' 




1 FGF2, AGGFl, CAVl, HIFIA, TNEgfi^^fi^UU^^HHI^^^I 




MIMAT0000077(hsa-miR-22-3p) 


GO term 


Genes 


Adj./>-value 


GO:0035583, Negative regulation of transforming growth factor 






MIMAT0000265(hsa-miR-204-5p) 


GO term 


Genes 


Adj.p-value 


GO:0035583, Negative regulation of transforming growth factor 
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Table S7. GOBO survival analysis results. Genes which was annotated with the specified 
GO term of proteins in the PIN would be used as input gene set for GOBO analysis. 
*p< 0.05, ** p< 0.01; *** p < 0.001. 



miRNA 



GO term 



/>-value 




GO:0015630, Microtubule cytoskeleton 

GO:0043067, Regulation of programmed cell death 
GO:0043069, Negative regulation of programmed cell death 
C^-mm-sAH Negative repulftinn r.f call H«»tli 



MIMAT0000076(hsa-miR-2 1 -5p) 




miR-22-3p) 



MIMAT0000089(hsa-miR-3 1 -5p) 



MAT0000250(hsa-miR- 1 39-5p) 



GO:0007173, Epidermal growth factor receptor signaling 
pathway 

GO:0043067, Regulation of programmed cell death 
GO:0035583, Negative regulation of transforming growth 
factor beta receptor signaling pathway by extracellular 
sequestering of TGFbeta 

GO:0003376, Sphingosine-1 -phosphate signaling pathway 
GO:0006917, Induction of apoptosis 

00:0007169, Transmembrane receptor protein tyrosine 
kinase signaling pathway 

GO:0007187, G-protein signaling, coupled to cyclic 
nucleotide second messenger 

GO:0007188, G-protein signaling, coupled to cAMP 
nucleotide second messenger 
GO:0010942, Positive regulation of cell death 
GO:0012502, Induction of programmed cell death 
GO:0019933, CAMP-mediated signaling 
GO:0019935, Cychc-nucleotide-mediated signaling 
GO:0043067, Regulation of programmed cell death 
GO:0043068, Positive regulation of programmed cell death 
GO:0048011, Nerve growth factor receptor signaling 
pathway 

GO:0035583, Negative regulation of transforming growth 
factor beta receptor signaling pathway by extracellular 
^sequestering of TGFbeta ^^^^^^^^^^^^^^^^^^m 



9.97 X lIHI 

0.268067 
0.0390439 * 
0.0390439 1 

0.721139 

0.266312 

0.940 



GO:0035583, Negative regulation of transforming growth 
MIMAT0000265(hsa-miR-204-5p) factor beta receptor signaling pathway by extracellular 

sequestering of TGFbeta 

GO:0002224, Toll-like receptor signaling pathway 
GO:0002755, MyD88-dependent toll-like receptor signaling 
pathway 

GO:0002756, MyD88-independent toll-like receptor 
signaling pathway 
GO:0006916, Anti-apoptosis 
GO:0006917, Induction of apoptosis 
iMAT0000423(hsa-miR-125b-5p) GO;0007169, Transmembrane receptor protein tyrosine 

kinase signaling pathway 

GO:0008284, Positive regulation of cell proliferation 
GO:0008285, Negative regulation of cell proliferation 
GO:0008629, Induction of apoptosis by intracellular signals 
GO:0010942, Positive regulation of cell death 




0.062202 
0.048584 * 

0.050408 
0.289466 

0.687572 

0.356228 
0.048584 * 
0.687572 
0.289466 
0.694486 
0.356228 

0.154543 
0.940727 J 



0.940727 

0.38C 
0.380928 




0.883 
0.883393 
0.073118 
0.972892 
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Table S7. Cont. 



miRNA 



GO term 



ju-value 




sa-miR-125b-5p) 



MIMAT0000437(hsa-miR-145-5p) 




IIMAT0000443(hsa-miR-125a-5p) 



GO:0015630, Microtubule cytoskeleton 

00:0030834, Regulation of actin filament depolymerization 
00:0030835, Negative regulation of actin filament 
depolymerization 

00:0030837, Negative regulation of actin filament 
polymerization 

00:0034130, Toll-like receptor 1 signaling pathway 
00:0034134, Toll-like receptor 2 signaling pathway 
00:0034138, Toll-like receptor 3 signaling pathway 
00:0034142, Toll-like receptor 4 signaling pathway 
00:0035666, TRIF-dependent toll-like receptor signaling 
pathway 

00:0043067, Regulation of programmed cell death 
00:0043068, Positive regulation of programmed cell death 
00:0043069, Negative regulation of programmed cell death 
00:0048011, Nerve growth factor receptor signaling 

pathway 

GO:0050678, Regulation of epithelial cell proliferation 

00:0050679, Positive regulation of epithelial cell 
proliferation 

00:0051693, Actin filament capping 

00:0060548, Negative regulation of cell death 

00:0007178, Transmembrane receptor protein 

serine/threonine kinase signaling pathway 

00:0030509, BMP signaling pathway 

00:0090092, Regulation of transmembrane receptor protein 

serine/threonine kinase signaling pathway 

00:0001936, Regulation of endothelial cell proliferation 

00:0002221, Pattern recognition receptor signaling pathway 

00:0002224, Toll-like receptor signaling pathway 

00:0002755, MyD88-dependent toll-like receptor signaling 

pathway 

00:0002756, MyD88-independent toll-like receptor 
signaling pathway 
00:0006916, Anti-apoptosis 
00:0006917, Induction of apoptosis 

00:0007169, Transmembrane receptor protein tyrosine 
kinase signaling pathway 

00:0008284, Positive regulation of cell proliferation 
00:0008285, Negative regulation of cell proliferation 
00:0008629, Induction of apoptosis by intracellular signals 
00:0010942, Positive regulation of cell death 
J3Q:0.Q12502..Induc.ticm.of programmed cell death 





0.654383 

0.380928 
0.380§, 
0.3 
0.3 

0.3 

0.643418 
0.972892 
0.492576 

0.171634 

0.002205 ** 

0.205483 

0.654383 
0.413665 

0.196953 

0.196953 

0.843529 

0.115146 
0.380928 
0.38092^ 

0.380928 




0.8? 

0.883393 
0.073118 
0.972892 J 



1^ 
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Table S7. Cont. 



miRNA 



GO term 



/j-value 



TMAT0000443(hsa-imR- 1 25a-5p) 




MIMAT0002820(hsa-miR-497-5p) 



GO:00 15629, Actin cytoskeleton 
GO:0015630, Microtubule cytoskeleton 

GO:0030834, Regulation of actin filament depolymerization 
GO:0030835, Negative regulation of actin filament 
depolymerization 

GO:0030837, Negative regulation of actin filament 
polymerization 

GO:0034130, Toll-like receptor 1 signaling pathway 
GO:0034134, Toll-like receptor 2 signaling pathw^ay 
00:0034138, Toll-like receptor 3 signaling pathway 
GO:0034142, Toll-like receptor 4 signaling pathway 
GO:0035666, TRIF-dependent toll-like receptor signaling 
pathway 

00:0043067, Regulation of programmed cell death 
00:0043068, Positive regulation of programmed cell death 
GO:0043069, Negative regulation of programmed cell death 
00:0048011, Nerve growth factor receptor signaling 
pathway 

GO:0050678, Regulation of epithelial cell proliferation 

00:0050679, Positive regulation of epithehal cell 
proliferation 

00:0050851, Antigen receptor-mediated signaling pathway 
00:0051693, Actin filament capping 
00:0060548, Negative regulation of cell death 
00:0001959, Regulation of cytokine-mediated signaling 
pathway 

00:0006917, Induction of apoptosis 

00:0007169, Transmembrane receptor protein tyrosine 
kinase signaling pathway 

00:0007173, Epidermal growth factor receptor signaling 
pathway 

00:0007178, Transmembrane receptor protein 
serine/threonine kinase signaling pathway 
GO:0007179, Transforming growth factor beta receptor 
signaUng pathway 

00:0008284, Positive regulation of cell proliferation 
00:0008285, Negative regulation of cell proliferation 
00:0008543, Fibroblast growth factor receptor signaling 
pathway 

00:0010942, Positive regulation of cell death 

00:0012502, Induction of programmed cell death 

00:0015629, Actin cytoskeleton 

00:0015630, Microtubule cytoskeleton 

00:0017015, Regulation of transforming growth factor beta 

receptor signaling pathway 



0.55 
0.028245 * 

0.654383 

0.6 




0.6 

0.38 
0.3? 

0.380928 
0.380928 

0.380928 

0.643418 
0.972892 
0.49257« 



0.1 

o.oo: 



0.205483 

0.103325 
0.6543J 




0.06699 
0.142401 
0.837635 

0.387447 

0.240167 

0.017876 * 

0.430255 
0.149994 

0.521978 

0.237692 
0.142401 
0.228804 
0.18331 

0.128773 
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Table S7. Cont. 



miRNA 



GO term 



^-value 




MIMAT0002856(hsa-miR-520d-3p) 



GO:0030509, BMP signaling pathway 
00:0030521, Androgen receptor signaling pathway 
GO:0032956, Regulation of actin cytoskeleton organization 
GO:0042058, Regulation of epidermal growth factor 
receptor signaling pathway 

00:0042059, Negative regulation of epidermal growth 

factor receptor signaling pathway 

00:0043067, Regulation of programmed cell death 

00:0043068, Positive regulation of programmed cell death 

00:0043069, Negative regulation of programmed cell death 

00:0048011, Nerve growth factor receptor signaling 

pathway 

00:0060548, Negative regulation of cell death 
00:0070302, Regulation of stress-activated protein kinase 
signaling cascade 

00:0090092, Regulation of transmemb rane receptor protein 
serine/threonine kinase signaling patl||[|^^^^^^^^^^^| 

00:0006917, Induction of apoptosis 

00:0007169, Transmembrane receptor protein tyrosine 
kinase signaling pathway 

00:0007173, Epidermal growth factor receptor signaling 
pathway 

00:0008284, Positive regulation of cell prohferation 
00:0008286, Insulin receptor signaling pathway 
00:0008543, Fibroblast growth factor receptor signaling 
pathway 

00:0008629, Induction of apoptosis by intracellular signals 
00:0010942, Positive regulation of cell death 
00:0012502, Induction of programmed cell death 
GO:0015630, Microtubule cytoskeleton 
00:0042058, Regulation of epidermal growth factor 
receptor signaling pathway 

00:0042059, Negative regulation of epidermal growth 

factor receptor signaling pathway 

00:0043067, Regulation of programmed cell death 

00:0043068, Positive regulation of programmed cell death 

00:0043069, Negative regulation of programmed cell death 

GO:0048011, Nerve growth factor receptor signaling 

pathway 

GO:0051988, Regulation of attachment of spindle 
microtubules to kinetochore 

00:0060548, Negative regulation of cell death 



0.3983" 
0.383811 
0.91762 

0.934041 

0.789493 

0.111544 
0.237691 
0.856891 

0.471981 

0.66743| 

0.561031 

^ 0.182314 
0.489781 
0.171689 

0.449696 

0.05916 
0.237933 

0.159318 

0.502822 
0.076906 
0.489781 
0.000292 *** 

0.762633 

0.826854 

0.740092 
0.076906 
0.067499 

0.014926 * 

7.48 X 10"* *** 
0.213035 
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Table S8. Pathophysiogical characteristics of miRNA array data used in ROC curve analysis. 



Sample Name 


ER 


PR 


HER 


TNM 


Stage 


Grade 


S621T1 


0 


0 


0 


pTlNOMO 


I 


2 


S434T1 


1 


0 


1 


T2N1M0 


IIB 


3 


S403T1 


0 


1 


0 


T2N1M0 


IIB 


2 


S459T1 


1 


0 


0 


T4N0M1 


IV 


1 


S455N1 


1 


0 


0 


T3N3M1 


IV 


3 


S545T1 


0 


0 


0 


pT2N0M0 


IIA 


3 


S173N1 


0 


0 


1 


T2N3M0 


IIIC 


3 


S363T1 


1 


1 


0 


T2N1M0 


IIB 


1 


S909T1 


1 


1 


1 


pT3N3aM0 


IIIC 


(Unknown) 


S645T1 


0 


1 


1 


pTlbNOMO 


I 


3 


S898T1 


1 


0 


0 


pT2N0(i-)M0 


IIA 


(Unknown) 


S201T1 


1 


0 


0 


TINIMO 


TT A 

IIA 


2 


S631T1 


1 


1 


1 


T2N3aMl 


IV 


2 


S303T1 


0 


0 


1 


T2N0M0 


IIA 


2 


S502T1 


0 


0 


0 


pT3N0M0 


IIB 


3 


S498N1 


1 


1 


1 


pTlcNlaMO 


IIA 


2 


S536T1 


1 


0 


1 


TlcNlmiMO 


IIA 


2 


S660T1 


0 


0 


1 


TINOMO 


I 


3 


S358N1 


0 


0 


1 


T2pN2M0 


IIIA 


2 


S665T1 


0 


0 


0 


T2N3M0 


IIIC 


3 


S475T1 


0 


0 


1 


pTlcNOMO 


I 


3 


S423T1 


1 


1 


0 


T2N3M0 


IIIC 


1 


S507T1 


0 


0 


0 


pTlcNlaMO 


IIA 


3 


S891T1 


1 


0 


0 


pT2NxM0 


IIA 


2 


S422T1 


1 


1 


0 


T2N0M0 


IIA 


1 


S961T1 


0 


0 


1 


T2N2aM0 


IIIA 


2 


S622T1 


0 


0 


0 


pT2N0M0 


IIA 


2 


S454T1 


0 


0 


0 


TlcNOMO 


I 


2 


S433T1 


1 


0 


0 


T2N0M0 


IIA 


1 


S673T1 


0 


0 


0 


pT2N0M0 


TT A 

IIA 


2 


S450T1 


0 


0 


1 


T2N1M0 


IIB 


3 


S430T1 


1 


1 


0 


T2N3M0 


IIIC 


2 


S437T1 


1 


0 


0 


T3N1M0 


IIIA 


3 


S574T1 


0 


0 


0 


T2N0M0 


IIA 


2 


S401T1 


1 


1 


0 


TINOMO 


I 


1 


CI /I T7XT 1 

S427N 1 


0 


0 


0 


T3N1M0 


IIIA 


2 


S894T1 


0 


0 


0 


pTlcNOMO 


I 


2 


S929T1 


1 


0 


0 


pT2N0M0 


IIA 


2 


S173T1 


0 


0 


1 


T2N3M0 


IIIC 


3 


S622N1 


0 


0 


0 


pT2N0M0 


IIA 


2 


S562T1 


1 


1 


1 


pTlaNOMO 


I 


3 


S602T1 


0 


0 


0 


pTlcNOMO 


I 


3 


S490T1 


1 


1 


1 


TT^XT/^A Jir\ 

pT2N0M0 


IIA 


2 


S677T1 


0 


0 


0 


pT2NlM0 


IIB 


3 


S881T1 


0 


0 


0 


pT2NlaM0 


IIB 


3 


S619T1 


0 


0 


0 


pT2N0M0 


IIA 


3 


S446T1 


0 


0 


0 


T2N2M0 


IIIA 


3 


S446T2 


0 


0 


0 


T2N2M0 


IIIA 


3 


S453T1 


1 


1 


1 


T2N1M0 


IIB 


2 


S562N1 


1 


1 


1 


pTlaNOMO 


I 


3 


S557T1 


0 


0 


0 


pT2N0M0 


IIA 


3 
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Table S8. Cont. 



Sample Name 




rR 


HER 


IJNM 


Stage 


Grade 


S594T1 


1 


1 


1 


pT2N3aM0 


IIIC 


3 


S582T1 


0 


0 


0 


pT2N0M0 


IIA 


3 


S358T1 


0 


0 


1 


T2pN2M0 


IIIA 


2 


S368T1 


0 


0 


1 


T3N3M0 


IIIC 


2 


S175T1 


1 


1 


0 


T3N3Mx 


IIIC 


3 


S357T1 


0 


0 


0 


T2pNlM0 


IIIC 


3 


S653T1 


0 


0 


0 


pT3N3aM0 


IIIC 


2 


S722T1 


1 


1 


1 


pTlNOMO 


I 


3 


S593T1 


0 


0 


0 


pTlNOMO 


I 


3 


S543T1 


1 


1 


1 


pT2NlM0 


IIB 


2 


S498T1 


1 


1 


1 


pTlcNlaMO 


IIA 


2 


S389T1 


1 


0 


0 


T2N1M0 


IIB 


3 


S614T1 


0 


1 


1 


TxNOMl 


IIIA 


(Unknown) 


S536N1 


1 


0 


1 


TlcNlmiMO 


IIA 


2 


S462T1 


1 


1 


0 


T3N3M0 


IIIC 


2 


S477T1 


0 


0 


0 


pTlcNOMO 


I 


3 


S917T1 


0 


0 


0 


T2N0M0 


IIA 


(Unknown) 


S213T1 


0 


0 


1 


TlcNlMO 


IIA 


3 


S628T1 


0 


0 


0 


TINOMO 


I 


2 


S291T1 


0 


0 


0 


T3pN2M0 


IIIA 


3 


S593N1 


0 


0 


0 


pTlNOMO 


I 


3 


S418T1 


1 


1 


0 


T2N3M0 


IIIC 


2 


S420T1 


1 


1 


0 


TINOMO 


I 


2 


S363N1 


1 


1 


0 


T2N1M0 


IIB 


1 


S629T1 


1 


0 


1 


pTlNOMO 


I 


2 


S586T1 


1 


1 


1 


pT2N3aM0 


IIIC 


2 


S415T1 


0 


1 


0 


T2N2M0 


IIIA 


2 


S439T1 


1 


1 


0 


T2N2M0 


IIIA 


2 


S941N1 


0 


0 


0 


pTlcNlaMO 


IIA 


3 


S893T1 


0 


0 


0 


pT2NlinicM0 


IIB 


3 


S380T1 


1 


0 


0 


T2N1mO 


IIB 


2 


f\r\ z' X T 1 

S906N1 


1 


1 


1 


pT3N3aM0 


IIIC 


2 


S918T1 


0 


0 


0 


pT2N0M0 


IIA 


3 


S400T1 


0 


0 


1 


T2N0M0 


IIA 


2 


S328T1 


0 


0 


0 


TlcpNOMO 


I 


3 


S367T1 


1 


0 


0 


TINIMO 


IIIA 


1 


S420N1 


1 


1 


0 


TINOMO 


I 


2 


S922T1 


0 


0 


0 


pTlcNO(i-)MO 


I 


3 


S896T1 


0 


0 


1 


pT2NlaM0 


IIB 


2 


S410T1 


0 


1 


1 


TINOMO 


I 


3 


S572T1 


0 


0 


1 


T4NlaM0 


IIIC 


3 


S448T1 


0 


0 


0 


TINOMO 


I 


3 


S207T1 


0 


0 


1 


T2N0M0 


IIA 


3 


S604T1 


0 


0 


1 


TT^X Tlx llf\ 

pT2NlM0 


IIB 


3 


S379T1 


0 


0 


1 


T2N2M0 


IIIA 


2 


S906T1 


1 


1 


1 


pT3N3aM0 


IIIC 


2 


S941T2 


0 


0 


0 


pTlcNlaMO 


IIA 


3 


S941T1 


0 


0 


0 


pTlcNlaMO 


IIA 


3 


S375T1 


1 


1 


0 


TINOMO 


I 


2 


S427T1 


0 


0 


0 


T3N1M0 


IIIA 


2 


S417T1 


0 


0 


0 


pT2N0M0 


IIA 


3 
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Table S8. Cont. 



Sample Name 


ER 


PR 


HER 


TNM 


Stage 


S180T1 


1 


0 


0 


T4NxM0 


IIIC 


S455T1 


1 


0 


0 


T3N3M1 


IV 


S887T1 




0 


0 


pT2N0M0 


IIA 


S445T1 


1 


0 


1 


T2N1M0 


IIB 


S469T1 


1 


1 


0 


T2N0M0 


IIA 


S469T2 


1 


1 


0 


T2N0M0 


IIA 


S483T1 






0 


T2N3M0 


IIIC 


S444T1 






0 


TlNxMO 


(Unknown) 


S909N1 






1 


pT3N3aM0 


IIIC 


S464T1 






0 


TINOMO 


I 


S894N1 


0 


0 


0 


pTlcNOMO 


I 


S698T1 


1 


1 


1 


pTlcNOMO 


I 


S452T1 


0 


0 


0 


TINOMO 


I 


S474T1 


(Unknown) 


(Unknown) 


(Unknown) 


(Unknown) 


(Unknown) 
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